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FOREWORD 


Simulating  the  interaction  of  terrain,  vehicle,  and  man  is 
essential  to  ing>roving  land  mobility  technology  — a goal  of  the 
U.S.  Army  Tank- Automotive  Command  (TACOM)  and  the  U.S.  Army  Corps 
of  Engineers  Waterways  Experiment  Station  (WES) . Iftider  TACOM  con- 
tract, Grumman  has  been  actively  supporting  this  effort  and  com- 
pleted the  development  of  a first  generation  mathematical  model 
for  driven  rigid  wheel  and  pneumatic  tire-soil  interaction. 

This  model,  conceived  as  an  alternate  submodel  in  the  terrain- 
vehicle-man  simulation,  called  the  AMC  Vehicle  Mobility  Model, 
allows  the  computation  of  tire  performance  in  both  purely  fric- 
tional or  cohesive  and  frictional-cohesive  soils. 

Ihis  first  generation  mathematical  model  is  based  on  the  con- 
cept of  "steady  state"  in  the  soil  that  assumes  a low  and  constant 
velocity  of  travel.  Observations  in  the  laboratory  and  in  the 
field  indicate  that  tire  performance  improves  with  travel  velocity 
under  certain  conditions. 

To  take  this  factor  into  account  in  mobility  evaluation  and 
take  advantage  of  it  in  off-road  vehicle  development  it  is  neces- 
sary to  anal3Tze  the  effects  of  soil  inertia  and  loading  rates  on 
tire-soil  interaction.  In  a follow-on  program  to  the  development 
of  a first  generation  pneumatic  tire-soil  model,  the  results  of 
which  are  presented  in  this  program,  the  model  was  expanded  to  in- 
clude soil  inertia  forces.  The  effects  of  the  loading  rate  were 
analyzed  by  adjusting  the  strength  properties  of  soil  for  the 
strain  rate  corresponding  to  the  travel  velocity  of  the  tire.  The 
expanded  model  was  also  shown  to  be  suitable  for  the  determination 
of  spring  constants  that  can  be  applied  in  vehicle  dynamic  models. 


For  the  analyses  of  cross  country  mobility  of  trailer-mounted 
weapons,  a towed  tire-soil  model  is  needed  that  predicts  the  motion 
resistance  and  sinkage  of  towed  pneumatic  tires.  Ihe  basic  con- 
cepts of  tire-soil  interaction  were  applied  to  the  case  of  towed 
tire  and  an.  appropriate  model,  described  in  the  second  part  of 
this  report,  has  been  developed. 


ABSTRACT 


The  role  of  soil  inertia  forces  and  the  effects  of  strain  rate 
in  tire-soil  interaction  are  analyzed.  Bie  differential  equations 
of  plasticity  for  soils  are  expanded  to  include  soil  inertia 
forces.  Two  methods  are  developed  for  the  determination  of  iner- 
tial accelerations  in  tire-soil  interaction.  In  the  first  method 
the  theory  of  velocity  fields  is  applied  for  the  determination  of 
accelerations  imparted  to  soil  particles  by  a tire  traveling  at 
constant  velocity.  An  iterative  scheme  is  developed  that  succes- 
sively updates  the  geometry  of  slip  line  fields  for  velocity  in- 
crements and  associated  inertia  forces.  In  the  second  method  soil 
particle  accelerations  are  computed  on  the  basis  of  particle  path 
geometry.  An  analytical  form  of  the  particle  path  is  assumed  and 
related  to  the  time  events  of  tire  travel.  Computations  by  both 
methods  indicate  that  the  solution  of  the  differential  equations 
of  plasticity  for  soil  becomes  multivalued  at  a relatively  low 
travel  velocity.  Up  to  this  velocity,  tire-soil  interaction  is  not 
significantly  affected  by  soil  inertia  forces.  Beyond  this  veloc- 
ity the  multivaluedness  of  solutions  of  the  differential  equations 
indicates  a condition  that  has  not  been  explored  yet  for  its  effect 
on  soil  behavior  and  cannot  be  treated  by  present  theories. 

The  effect  of  strain  rate  of  soil  strength  properties  is  ana- 
lyzed and  recognized  as  the  major  contributor  to  the  improvement 
of  tire  performance  with  speed. 

A mathematical  model  for  towed  pneumatic  tires  and  soil  is  de- 
veloped for  the  simulation  of  towed  tire-soil  interaction.  Towed 
force  coefficient  predictions  are  compared  with  experimental  re- 
sults . 
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I.  SCOPE  OF  WORK 


The  scope  o>:  work  as  described  in  the  RFP  work  statements  for 
Contracts  DAAE07 -74-Q-0004  and  C-0204  was  twofold:  investigate 

the  effect  of  spe>.d  on  tire-soil  interaction  and  develop  a towed 
pneumatic  tire-soil  model.  Within  this  general  scope  emphasis  was 
placed  on  the  following  items: 

• Expansion  of  plasticity  theory  for  soils  to 
include  soil  inertia  forces 

• Development  of  methods  for  the  determination 
of  inertia  accelerations  of  soil  particles 
caused  by  tire  travel 

• Analyses  of  soil  inertia  and  strain  rate  ef- 
fects on  tire-soil  interaction 

• Development  of  a towed  tire -soil  model  with 
consideration  of  interface  stress  distributions 
resulting  in  zero  torque  condition 

• Verification  of  the  towed  tire-soil  model  by 
comparisons  of  predicted  towed  force  coeffi- 
cients with  available  experimental  data. 


1 


II.  EFFECT  OF  SPEED  ON  TIRE-SOIL  INTERACTION 


At  present,  tire-soil  interaction  theories  are  based  on  a 
quasi-static  model  that  assumes  a steady  state  of  the  tire-soil 
system.  Ihe  tire  is  assumed  to  travel  at  a sufficiently  low 
velocity  so  that  velocity  effects  are  negligible.  Yet  experiments 
performed  under  controlled  conditions  at  WES  (Ref.  1)  show  that 
even  in  the  low  velocity  range  of  off-road  travel,  velocity  ef- 
fects are  far  from  negligible.  In  the  higher  velocity  range  of 
aircraft  landing  gear-soil  interaction  velocity  is  one  of  the  main 
controlling  factors  that  affects  the  drag  load  as  experiments  indi- 
cate (Refs . 2 and  3) . 

Velocity  affects  tire-soil  interaction  primarily  in  the  fol- 
lowing two  ways: 

• Soil  inertia  forces  are  generated  in  the  soil 
during  the  passage  of  a tire.  The  magnitude  of 
these  inertia  forces  depends  on  the  velocity  of 
the  tire  and  is  approximately  proportionate  with 
the  square  of  the  velocity  of  the  tire. 

• The  strength  of  the  soil  that  controls  the  inter- 
face stresses  is  affected  by  the  rate  of  loading 
that  is  directly  proportionate  to  the  velocity 

of  travel. 

Secondary  effects  of  velocity  on  tire-soil  interaction  may  in- 
clude the  increase  of  tire  stiffness  with  an  increase  of  the  angu- 
lar velocity  of  tire,  the  effects  of  velocity  on  interface  friction 
development,  etc.  No  study  of  these  effects  was  made  in  the  pres- 
ent research  project. 
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III.  EFFECT  OF  SOIL  INERTIA  FORCES 
ON  TIRE-SOIL  INTERACTION 


The  concept  of  tire-soil  interaction  presented  in  Ref.  4 for 
driven  tires  and  Section  VIII  of  this  report  for  towed  tires 
assumes  that  the  soil  is  in  the  plastic  state  of  equilibrium,  char- 
acterized by  slip  line  fields,  whenever  the  normal  stress  corre- 
sponding to  this  plastic  state  of  stresses  is  less  than  the  limit- 
ing pressure  characteristic  of  the  tire.  The  geometry  of  the 
zones  of  plastic  equilibrium  and  the  stress  states  in  these  zones 
are  defined  by  solutions  of  the  differential  equations  of  plas- 
ticity for  soils.  These  differential  equations  are  derived  from 
the  combination  of  the  differential  equations  of  static  equilibrium 
with  the  Mohr-Coulomb  yield  criterion. 

The  passage  of  a tire,  even  if  it  travels  at  a constant  veloc- 
ity, displaces  soil  particles  within  the  affected  soil  mass.  The 
end  result  of  soil  particle  displacement  is  the  formation  of  a rut 
behind  the  tire.  Thus,  during  the  passage  of  the  tire  a soil 
particle  moves  from  its  original  position  at  rest  to  a new  one 
where  it  comes  to  rest  again.  This  particle  motion  involves  ac- 
celeration and  deceleration  of  the  particle  and,  as  a consequence, 
inertia  forces  are  generated  in  the  soil.  To  consider  these  inertia 
forces  in  the  differential  equations  of  plasticity  for  soil,  the 
Mohr-Coulomb  yield  criterion  has  to  be  combined  with  the  differen- 
tial equations  of  motion  instead  of  those  of  static  equilibrium. 

This  combination,  after  the  same  manipulation  of  the  equations  is 
performed  as  in  the  static  case,  yields  the  following  differential 
equations: 


dz  = dx  tan(0  ^ (j.) 


(1) 


da^  Za  tan  <pd0  = — 

S 


^a^^  (a^  + g)tan  <p)dx+(a^+g  ± a^  tan  <p)dz| 


s Preceding  page  blank 


In  these  equations  x and  z refer  to  a coordinate  system 

that  moves  with  the  tire,  whereas  a and  a denote  soil  parti- 

X z 

cle  accelerations  that  refer  to  a fixed  coordinate  system.  Ihe 
motion  of  soil  particles  as  well  as  their  accelerations  are  closely 
related  to  the  geometry  of  slip  line  fields.  Since  this  latter 
itself  depends  on  the  particle  accelerations,  as  Eq.  (1)  indicates, 
therefore,  there  is  an  interdependence  between  the  accelerations 
and  other  variables  in  Eq.  (1).  Even  if  a functional  relationship 
could  be  established  among  these  variables,  the  direct  numerical 
solution  of  Eq.  (1)  would  be  a formidable  task  because  of  the  in- 
herent nonlinearities  in  these  relationships.  In  an  attempt  to 
develop  a method  that  yields  an  acceptable  approximate  solution  of 
Eq.  (1),  two  approaches  have  been  followed  and  investigated  in  de- 
tail. In  the  first,  essentially  theoretical  approach  the  theory 
of  velocity  fields  associated  with  the  slip  line  fields  is  applied 
to  the  determination  of  particle  accelerations.  Ihe  second  ap- 
proach is  based  on  an  analytical  simulation  of  observed  particle 
path  geometries  and  is,  therefore,  a semi-empirical  one.  A de- 
tailed discussion  of  these  approaches  together  with  an  evaluation 
of  the  numerical  computational  procedures  developed  for  their  ap- 
plication to  the  tire-soil  interaction  problem  is  given  in  the  fol- 
lowing sections. 
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IV.  DETERMINATION  OF  SOIL  INERTIA  FORCES 
BY  THE  THEORY  OF  VELOCITY  FIELDS 


Theoretical  Background 

Originally,  the  theory  of  velocity  fields  was  developed  in 
connection  with  the  application  of  the  theory  of  plasticity  to 
metal  forming  processes.  For  metals  the  angle  of  internal  fric- 
tion is  zero  and  the  characteristics  obtained  from  the  solutions 
of  the  differential  equations  for  both  the  velocities  and  the 
stresses  coincide.  In  the  plastic  state,  frictionless  materials 
actually  slip  along  the  characteristic  line;  hence  the  term  "slip 
line  field." 

The  theory  of  velocity  fields  was  later  extended  for  materials 
exhibiting  friction,  such  as  soils.  The  purpose  of  this  extension 
was,  however,  not  so  much  to  determine  the  flow  of  the  material, 
as  in  the  case  of  metal  forming  processes,  but  to  enable  research- 
ers to  apply  the  limit  theorems  of  plasticity  to  various  problems 
of  the  plastic  equilibrium  in  soils.  The  limit  theorems  of  plas- 
tic equilibrium  state  that  a statically  admissible  stress  field, 
derived  from  the  differential  equations  of  the  plasticity  theory, 
yields  a lower  bound  for  the  collapse  load,  while  collapse  loads 
computed  for  various  kinematically  admissible  velocity  fields 
yield  upper  bound  solutions.  The  true  solution  is  when  the  lower 
and  upper  bounds  coincide. 

The  significance  of  the  limit  theorems  of  plasticity  in  engi- 
neering applications  is  that  the  collapse  load  (or  failure  stresses) 
may  be  higher  than  that  computed  from  stress  fields  if  the  geometry 
of  the  stress  fields  is  kinematically  inadmissible. 

Derivation  of  the  differential  equations  of  velocity  fields 
is  based  on  the  assumption  that  the  strain  rate  is  proportional  to 


7 


the  time  rate  of  stress  (Refs.  5 through  12).  A system  of  partial 
differential  equations  similar  to  that  of  the  stress  equations  can 
then  be  derived.  The  solutions  of  this  set  of  hyperbolic  differen- 
tial equations  are  the  characteristics  of  the  velocity  field.  Gen- 
erally, the  kinematic  boundary  conditions  are  insufficient  to  define 
the  velocity  field  and,  therefore,  the  usual  procedure  is  to  ana- 
lyze the  geometries  of  stress  fields  for  kinematic  admissibility. 

If,  however,  for  a frictional  material  the  characteristics  of  the 
velocity  field  are  assumed  to  coincide  with  those  of  the  stress 
field,  the  strain  rates  associated  with  the  velocity  fields  require 
a rate  of  dilation  of  the  soil  that  is  unrealistic.  Various  theo- 
ries have  been  proposed  to  deal  with  this  problem  but  none  of  them 
is  completely  satisfactory.  It  is  outside  the  scope  of  this  report 
to  discuss  these  controversial  issues  in  detail;  regarding  the  use 
of  velocity  fields  for  the  determination  of  soil  inertia  forces  in 
tire-soil  interaction  problems  the  following  considerations  apply: 

• If  the  stress  boundary  conditions  are  such  that 
they  uniquely  determine  the  stress  field  then 
the  problem  is  statically  determinate,  the  char- 
acteristics of  the  velocities  coincide  with  those 
of  the  stresses.  The  solution  is  the  true  solu- 
tion since  upper  and  lower  bounds  are  identical 
(Ref.  11). 

• In  tire-soil  interaction  problems  boundary  condi- 
tions completely  define  the  stress  field,  there- 
fore, the  solution  is  statically  determinate  and 
it  is  the  true  solution.  However,  the  assumption 
of  a uniform  development  of  interface  friction 
along  the  contact  area,  although  to  some  extent 
supported  by  experiments,  is  an  arbitrary  one. 
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Since  interface  friction  depends  on  the  kinematic 
boundary  conditions  at  the  interface,  the  solution 
of  Eq.  (1)  is  the  true  solution  only,  if  the  as- 
sumed interface  friction  distribution  is  compati- 
ble with  the  kinematic  boundary  conditions.  The 
establishment  of  the  kinematic  boundary  conditions 
at  the  interface  is  also  necessary  for  the  compu- 
tation of  velocities  and  is  discussed  in  detail 
under  that  heading. 

• In  the  following  analyses  it  is  assumed  that  the 
characteristics  of  stress  and  velocity  fields 
coincide.  Although  this  assumption  is  associated 
with  an  unrealistic  rate  of  dilation,  in  the  case 
of  tire-soil  interaction  problems  the  velocity 
field  is  an  instantaneous  one,  and  the  dilation 
prevails  only  for  an  infinitesimal  time.  Also, 
the  tire-soil  models  used  for  the  calculation  of 
tire  performance  assume  a steady  state  in  soil 
where  the  soil  is  already  compacted  by  the  tire 
action  to  a significant  degree.  Therefore,  the 
soil  is  likely  to  dilate  under  the  yield  stresses. 

In  sunmary,  flow  fields  in  tire-soil  interaction  are  main- 
tained for  an  infinitesimal  time,  therefore,  the  rate  of  dilation 
is  not  a critical  consideration  and  the  characteristics  of  stress 
and  velocity  fields  may  be  assumed  to  coincide.  If  the  geometry 
of  the  characteristics  is  computed  from  the  equations  for  the 
stress  field  then  the  velocities  can  be  computed  from  the  follow- 
ing equations  (Ref.  10). 
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(2) 


dv  + (v  tan  cp  + v sec  <p)d0  * 0 

dv^  - (v^  sec  qp  + tan  cp)d0^  = 0 

The  above  equations  define  the  variation  of  velocities  along 
characteristic  lines.  The  velocity  vectors  v and  v^  are  com- 
ponents (projections)  of  the  velocity  vector  v in  the  direction 
of  the  "j"  and  "i"  lines,  respectively.  It  is  noted  that  the 
vectorial  sum  of  v and  v^  is  not  v (Fig.  1). 

The  V and  v components  of  the  velocity  vector,  v,  are 

X Z >s 

ct  6 

related  to  v and  v^  as  follows 

ct  . 

V = V sin  e - V cos  e 

X z 

(3) 

& . & 6 

V = V sin  0 - V cos  0 

X z 

Since  it  is  assumed  that  the  velocity  and  stress  characteris- 
tics coincide,  the  angles  0^  and  0^  may  be  determined  from  the 
stress  field  relation 

= e + (f  - D 

(4) 

For  the  numerical  computation  of  the  velocities  the  geometry 
of  the  characteristic  lines  as  well  as  the  0 values  are  computed 
for  the  stress  field  and  the  velocities  from  the  following  differ- 
ence equations; 

and  P are  superscripts  designating  slip  lines 
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(A  - B • C)/(1.0  - B • D) 


where 


a 

V.  . 


V?  . 

3-0 


« « « 

C - D • V.  . 


A = 


1 + i tan  <P(9.^.  - 


i sec  - 9^41^,) 

1 + i tan  ?(9  - 9 


c = 


Ij-H  -*•  *n,j+l  <^(^1,3  - ^l,J+l) 

1 - i tan  q>(9.^^  - 9.^ 


D » 


i 'l’<^l,j  - ^l,J+l) 
1 - i tan  <p(9^^j  - 


(5) 


The  difference  equations  (5)  indicate  that  for  given  stress 
field  characteristics  the  velocity  for  a grid  point  i,j  can  be 
computed  if  the  velocities  at  two  adjacent  points  (i+l,j  and 
i,j+l)  are  known.  The  kinematic  boundary  conditions,  discussed 
under  the  next  heading,  provide  the  initial  values  of  velocities 
necessary  to  start  the  numerical  computations. 


Kinematic  Boundary  Conditions 

In  a slip  line  field  determined  by  the  differential  equations 
of  plasticity  for  soils  by  methods  discussed  in  detail  in  Refs.  13 
through  18  there  are  three  boundaries  where  kinematic  boundary  con- 
ditions are  to  be  examined  (Fig.  2).  At  the  stress  free  surface 
AB  it  is  assumed  that  there  are  no  kinematic  constraints  and 
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points  at  the  surface  may  be  displaced  freely.  The  boundary  BCDE 
is  the  outermost  slip  line  beyond  which  the  soil  is  not  in  the 
state  of  failure.  Implicit  in  the  use  of  the  Mohr-Coulomb  yield 
criterion  is  that  soil  deformations  prior  to  yield  are  disregarded, 
therefore,  in  this  concept  the  soil  beyond  the  outer^st  slip  line 
acts  as  a rigid  body.  Consequently,  velocities  across  this  bound- 
ary must  be  zero,  resulting  in  the  kinematic  boundary  condition 
= 0 for  this  line.  For  this  boundary  condition  Eq.  (~)  is 
integrable  and  yields 


a 


V 


. ct  a. 
tan  <p(0  -0  ) 

e 


(6) 


Note  that  for  the  stress  computations  the  state  of  the  soil 
outside  the  slip  line  field  is  immaterial.  The  velocity  field, 
however,  is  profoundly  affected  by  the  assumption  that  the  soil 
adjacent  to  the  outermost  slip  line  is  assumed  to  be  rigid.  Ex- 
perimental information  on  soil  displacements  along  this  boundary 
is  scarce.  It  appears  that  while  the  soil  mass  outside  this 
boundary  undoubtedly  undergoes  some  deformation  there  is  an  abrupt 
change  in  the  displacements  in  a narrow  band  along  the  boundary. 
Thus,  the  theoretical  boundary  condition  corresponding  to  the 
rigid-plastic  material  idealization  may  not  be  too  far  off  from 
reality. 

As  stated  before,  the  critical  kinematic  boundary  conditions 
are  those  at  the  tire-soil  interface.  The  distribution  of  interface 
friction  along  this  boundary  must  be  consistent  with  the  kinematic 
conditions  there  if  the  slip  line  field  solution  is  to  be  the  true 
solution  for  the  yield  stresses.  As  a result  of  detailed  studies 
of  the  kinematic  conditions  at  the  interface,  new  concepts  of  slip 
and  interface  friction  were  developed.  These  are  not  only  useful 
for  the  definition  of  boundary  conditions  for  velocity  fields  but 
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also  offer  new  insight  in  the  not  very  well  understood  area  of 
si ip -interface  friction  relationships.  Hie  results  of  these  kine- 
matic analyses  and  the  reasoning  that  led  to  the  development  of 
these  new  concepts  regarding  interface  friction  are  presented  in 
the  following  discussions. 

The  kinematics  of  a point  attached  to  the  surface  of  the  tire 
are  easily  determined  if  the  geometry  of  the  tire  centerline  and 
the  angular  and  translational  velocities  of  the  tire  are  known. 

Hie  trajectory  of  such  a point,  in  a fixed  coordinate  system  is 
shown  in  Fig.  3 for  a complete  revolution  of  a slipping  tire.  Hie 


Fig,  3.  Trajectory  of  a Point  at  the 

Circumference  of  a Deformable  Tire 

point  describes  a prolate  cycloid  when  its  position  is  in  the  un- 
deflected part  of  the  tire;  otherwise  the  trajectory  is  an  irregu- 
lar curve  corresponding  to  the  assumed  tire  deflection.  The  hori- 
zontal and  vertical  components  of  the  velocity  of  the  point  may 
also  be  easily  determined  since  both  the  angular  and  translational 
velocities  are  assumed  to  be  constants. 

For  the  determination  of  the  velocity  boundary  conditions 
that  apply  to  the  slip  line  field  in  the  soil  the  velocities  im- 
parted to  a soil  particle  at  the  interface  are  needed.  These  are 
related  to  the  velocities  of  a point  attached  to  the  tire  surface 
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but  are  by  no  means  identical.  A fundamental  aspect  of  the  analy- 
sis of  the  kinematic  conditions  at  the  interface  is  that  the  inter- 
face must  be  considered  as  consisting  of  two  faces;  one  being  the 
surface  of  the  tire,  the  other  that  of  the  soil.  It  is  assumed 
that  these  two  faces  of  the  interface  may  slide  on,  but  may  not 
separate  from  each  other.  This  latter  assumption  leads  to  the 
condition  that  the  normal  components  of  the  velocities  of  the  two 
faces  be  the  same  at  every  point.  This  condition,  however,  is  by 
itself  insufficient  for  the  determination  of  soil  particle  veloci- 
ties at  the  interface.  Hie  additional  assumption  is  made  that  the 
direction  of  the  velocity  vector  of  a soil  particle  coincides  with 
the  direction  of  the  major  principal  stress  in  the  slip  line  field. 
For  homogeneous  isotropic  soils  with  linear  stress-strain  relations 
this  assumption  follows  from  the  theorems  of  continuum  mechanics; 
for  soils  not  meeting  the  above  criteria  strictly  the  assumption 
is  only  approximately  valid.  With  these  assumptions  the  velocities 
of  soil  particles  at  the  interface  (that  constitute  the  velocity 
boundary  conditions)  can  be  determined  as  follows. 

Ihe  velocity  vector  for  a point  at  the  tire  surface  may  be  de- 
termined by  vector ially  adding  the  forward  velocity  vector  (v^ 
and  the  tangential  velocity  vector  (v^)  (Fig.  4a).  The  tangen- 
tial velocity  vector  is  the  component  of  vector  v (peripheral 
velocity)  in  the  direction  of  the  deflected  surface  of  the  tire 


V 

P 


cor  = v^/  (1  - s) 


V 


t 


V cos  (a 
P 


a') 


(7) 


The  horizontal  and  vertical  components  of  the  resultant  veloc- 
ity vector  are 
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editions  at  the  Tire-Soil  Interface . 
For  Soil 
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V-  - V cos  (a  - a')  cos  a* 
f p ' 


V =*  V cos  (a  - a') sin  a' 
z p 

V « V sin  + V cos  a'  * v-  sin  a' 

n X z f 


(8) 


This  last  equation  expresses  the  fact  that  the  con^onent  of 
the  velocity  vector  nonnal  to  the  deflected  surface  of  the  tire 
depends  only  on  the  forward  velocity  and  is  independent  of  slip. 

For  a soil  particle  at  the  interface  the  nonnal  component  of 
the  velocity  vector  can  be  computed  as  (Fig.  4b) 

t 

S S./.S  / ./ 

V - V sin  a + V cos  a = v-  sin  a 
n X z f 

(9) 


= tan  6' 


From  Eq.  (8)  the  horizontal  and  vertical  components  of  the 
velocity  vector  can  be  expressed  as 


s 

V 

X 


sin 

'^f  sin  + cotan  9*  cos  a' 


s s . 

V = V cotan  Q 
Z X 


(10) 


The  velocity  vector  of  a soil  particle  at  the  interface  can 
also  be  thought  of  as  being  composed  of  the  forward  velocity  vector 
and  a tangential  velocity  vector.  This  latter  can  be  computed  from 
the  fol'ioiring  equation: 


s 1 

'^t  ^f  a'  tan  9'  + cos  a' 


(11) 


Superscript  s 


designates  soil  particle  velocities 
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The  tangential  velocity  conqjuted  by  Eq.  (11)  for  a soil  parti- 
cle is  in  the  case  of  driven  tires  generally  less  than  the  tangen- 
tial velocity  computed  by  Eq.  (7)  for  a point  at  the  surface  of 
the  tire  indicating  that  a relative  displacement  occurs  between 
the  soil  and  the  tire.  From  the  value  of  the  tangential  velocity 
of  a soil  particle  a hypothetical  slip  value  may  be  computed  as 

s'  = 1 - cos  (a  - a')  ♦ (sin  a*  tan  0*  + cos  a')  (12) 

It  is  interesting  to  note  that  this  h3q>othetical  slip  value 
is  independent  of  the  translational  velocity  of  the  tire,  at  least 
for  the  assumed  steady  state  concept  of  tire-soil  interaction. 

At  this  point  it  is  worthwhile  to  consider  the  meaning  of  the 
h3rpothetical  slip  value  determined  by  Eq.  (12).  If  the  tire  ac- 
tually turned  with  an  angular  velocity  that  resulted  in  a slip 
equal  to  s ' then  there  would  be  no  differential  displacement  be- 
tween the  surface  of  the  tire  and  the  soil.  Therefore,  the  hypo- 
thetical slip  s'  may  be  termed  as  "contact  slip."  On  the  other 
hand,  the  difference  between  s and  s'  is  the  result  of  dif- 
ferential displacement  between  tire  surface  and  soil.  The  term 
"spin"  conveys  the  idea  of  a tire  rotating  relative  to  a nonyield- 
ing surface  and  is  also  expressive  of  the  relative  motion  between 
a tire  and  a yielding  surface,  therefore  it  will  be  used  herein  to 
designate  the  quantity  s - s'. 

To  further  elucidate  the  meaning  and  significance  of  these 
terms  a track  element  may  be  considered  (Fig.  5a).  If  a track  ele- 
ment displaces  the  soil  by  AL  but  stays  in  contact  with  the  soil 
while  the  vehicle  travels  a distance  L then  the  "contact  slip"  is 

s'  = (l  - ^ (13) 
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Pig.  5.  Contact  Slip  (a)  and  Total 
Slip  (b)  of  a Track  Element 

However,  if  the  track  element  slides  over  the  soil  and  at  the 
end  of  vehicle  travel  L is  displaced  by  AX  relative  to  its 
original  position  (Fig.  5b)  then  the  total  slip  is 

s = (14) 

In  the  case  of  a track  element  it  is  easy  to  visualize  the 
meaning  of  the  two  slip  values  s and  s'.  In  the  case  of  tire- 
soil  interaction  the  assumption  of  a steady  state  in  the  soil 
leads  to  the  concept  that  the  slip  line  fields  move  with  the  tire 
while  the  soil  particles  are  displaced  from  their  original  position 
with  respect  to  a coordinate  system  that  is  fixed  to  the  ground. 
This  concept  makes  it  somewhat  difficult  to  visualize  the  relative 
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displacement  of  the  tire  surface  to  the  soil.  However,  the  meaning 
of  the  terms  introduced  by  Eq.  (12)  is  essentially  the  same  as  ex- 
plained for  a track  element. 

In  the  foregoing  discussions  velocities  of  a point  at  the  sur- 
face of  tire  and  soil  were  defined  and  compared.  On  the  other  hand, 
the  conventional  definition  of  slip  relates  slip  to  actual  and 
theoretical  distances  traveled  by  the  tire.  Obviously,  the  con- 
ventional slip  value  determines  an  overall  slip  for  the  tire. 

Janosi  (Ref.  19)  and  others  have  pointed  out  that  the  slip  for  an 
infinitesimal  element  of  a rigid  wheel  may  vary  along  the  interface. 
In  the  case  of  pneumatic  tires  infinitesimal  elements  of  the  sur- 
face of  both  tire  and  soil  may  deform  and  exhibit  different  dis- 
placements relative  to  each  other.  As  a result,  a variation  of 
slip  along  the  interface  may  easily  develop. 

The  kinematic  boundary  conditions  and  Eq.  (12)  derived  from 
these  conditions  allow  the  computation  of  contact  slip  for  every 
point  at  the  interface  once  the  tire  centerline  geometry  and  the 
associated  slip  line  fields  have  been  determined.  For  the  determi- 
nation of  the  slip  line  fields  in  the  tire-soil  model  the  angle  of 
interface  friction  is  assumed,  or  computed  from  the  modified 
Janosi-Hanamoto  equation 


tan  6 = 


tan  6 


max 


(1 


-(s+s^)/k^ 


(15) 


for  given  slip  and  slip-shear  parameters.  Equation  (15)  refers  to 
the  conventional  value  of  slip.  In  the  model,  slip  and  interface 
friction  angle  are  assumed  to  be  constant  along  the  interface. 

For  these  conditions  the  variation  of  the  contact  slip  along  the 
interface  in  the  tire-soil  model  was  determined  for  several  loading 
conditions,  tire  t3q)es,  and  conventional  slip  values.  It  was  found 


20 


that  the  variation  of  the  contact  slip  along  the  interface  was 
minimal  and  its  value  less  than  the  conventional  value  of  slip. 

T3TJical  variations  of  the  contact  slip  along  the  interface 
are  shown  in  Fig.  6 as  a function  of  the  central  angle.  The  con- 
tact slip  values  were  found  to  vary  little  along  the  interface  in 
the  cases  investigated.  Largest  deviations  from  the  average  value 
were  found  near  the  exit  angle  where  the  tire  centerline  geometry 
shows  an  upward  curvature.  Kinematically  this  portion  of  the 
interface  is  the  most  problematic  since  the  kinematic  boundary 
conditions  as  set  forth  in  the  foregoing  paragraphs  require  that 
the  velocity  of  the  soil  particle  be  upward  directed  so  that  no 
separation  of  soil  and  tire  surface  occurs. 

The  concept  of  contact  slip  has  important  implications  for 
both  the  theory  of  tire-soil  interaction  and  the  practical  aspects 
of  off-road  vehicle  engineering.  The  contact  slip  represents  the 
absolute  minimum  of  slip  at  which  a certain  tire  performance  can 
be  realized.  In  the  tire-soil  model  the  tire  centerline  geometry 
as  well  as  the  directions  of  the  principal  stresses  along  the  in- 
terface are  computed  and  from  this  information  contact  slip  values 
for  every  nodal  point  at  the  interface  can  be  determined.  Since 
the  variation  for  the  contact  slip  values  along  the  interface  is 
minimal,  an  average  value  of  the  contact  slip  may  be  computed  and 
relationships  between  the  development  of  interface  friction  angle 
and  the  average  value  of  the  contact  slip  may  be  established  for 
various  conditions.  Thus  a theoretical  basis  exists  for  analysis 
of  the  relationship  between  interface  friction  angle  and  contact 
slip.  This  relationship  may  be  compared  with  Eq.  (15)  that  repre- 
sents the  relation  with  the  conventional  slip  value. 

The  concept  of  contact  slip  and  the  method  of  analysis  that 
allows  its  determination  for  various  soil  and  loading  conditions 
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opens  new  avenues  of  research  that  could  lead  to  economics  in  off- 
road vehicle  engineering.  Obviously,  spin,  the  differential  dis- 
placement between  tire  and  soil  surface,  is  undesirable  since  it 
is  the  primary  cause  of  tire  wear  and  results  in  a waste  of  energy. 
For  understanding  of  the  phenomena  that  causes  spin  a systematic 
analysis  of  the  contact  slip  and  its  relation  to  the  conventional 
slip,  soil,  and  loading  conditions, preferably  coupled  with  ex- 
perim<=^ntal  investigations,  would  be  needed.  These  would  include 
measurements  of  soil  particle  displacements  at  the  surface  as  well 
as  torque  rather  than  slip-controlled  perfomiance  tests. 

Computation  of  Inertial  Accelerations  from  Velocity  Fields 

The  concept  of  tire-soil  interaction  for  which  the  mathemati- 
cal model  was  developed  (Ref.  13)  assumes  that  steady  state  condi- 
tions exist  in  the  soil.  This  assumption  is  equivalent  to  stating 
that  the  slip  line  fields  travel  with  the  tire  with  the  same 
velocity.  For  thes<i  conditions  the  computation  of  accelerations 
becomes  very  simple  if  the  velocity  fields  are  known.  In  an  in- 
finitesimal At  time  increment  the  velocity  of  a soil  particle  at 
any  point  x,z  changes  to  the  velocity  at  x - Ax,z  where 
Ax  = v^t.  The  velocity  field  computation  yields  the  velocity  com- 
ponents for  every  nodal  point  in  the  slip  line  field.  The  numeri- 
cal computation  of  the  accelerations  for  these  nodal  points  re- 
quires the  determination  of  the  velocities  for  a Ax  change.  This 
is  accomplished  by  assuming  a linear  velocity  change  between  the 
nodal  point  in  question  and  two  adjacent  nodal  points,  one  along 
the  "i"  and  one  along  the  "j"  line.  At  the  boundaries  of  the 
slip  line  field,  velocities  at  these  points  are  not  always  avail- 
able and  other  combinations  of  adjacent  nodal  points  had  to  be 
selected  for  the  determination  of  accelerations. 
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Computational  Scheme  for  the  Consideration  c£ 

Soil  Inertia  Forces  in  Tire-Soil  Interaction 

For  analysis  of  the  effect  of  soil  inertia  forces,  computed 
on  the  basis  of  velocity  fields,  a computer  program  was  prepared. 
This  computer  program  computes  the  interface  stresses  for  a given 
entry,  rear  and  separation  angles  for  various  translational  veloc- 
ities of  the  tire.  The  soil  inertia  forces  considered  in  the 
program  correspond  to  a constant  translational  velocity  of  the 
tire. 

The  program  represents  a numerical  solution  method  for  the 
differential  equations  (1) . It  computes  the  soil  inertia  forces 
by  successive  iteration.  In  this  iteration  the  first  step  is  to 
determine  the  slip  line  field  geometry  without  soil  inertia  forces, 
i.e.,  for  v^  = 0.  Then  on  the  basis  of  this  slip  line  field 
geometry  the  velocity  fields  and  accelerations  are  determined  for 
each  nodal  point  for  a selected  velocity  increment.  It  is  assumed 
that  the  accelerations  computed  for  each  nodal  point  are  valid 
when  a new  slip  line  field  geometry  is  computed  on  the  basis  of 
these  accelerations.  This  iteration  scheme  in  the  computer  pro- 
gram is  included  in  an  interactive  manner  that  allows  study  of  the 
effect  of  the  selection  of  velocity  increments  as  well  as  the 
effect  of  repeated  iterative  acceleration  computations  on  the  slip 
line  field  geometry.  A detailed  flow  chart  for  the  program  is 
given  in  Appendix  A. 

Problems  Encountered  with  Development  of 
the  Computer  Program 

Numerous  minor  problems  were  encountered  with  development  of 
the  program  because  initially  the  same  iteration  and  interpolation 
schemes  were  used  in  the  program  as  in  the  case  without  inertia 
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forces.  TSie  consideration  of  inertia  forces  in  Eq.  (1)  resulted 
in  a more  complex  situation  where  the  behavior  of  soil  was  often 
different  from  that  anticipated  on  the  basis  of  previous  experi- 
ence and,  therefore,  iteration  and  interpolation  schemes  had  to  be 
modified.  Other  problems  that  have  broader  significance  are  dis- 
cussed below. 

The  starting  point  for  the  development  of  this  computer  pro- 
gram was  the  one  prepared  for  the  simulation  of  tire-soil  interac- 
tion presented  in  Ref.  13.  In  that  program  an  ixj  --  48  x 16 
grid  is  set  up  for  the  computation  of  a slip  line  field,  but  only 
as  many  j lines  are  computed  as  are  necessary  to  end  up  with  the 
boundary  j line  at  a prescribed  location  at  the  interface.  The 
boundary  j line  is  obtained  by  interpolation  between  evenly 
spaced  j lines.  This  technique  was  developed  to  reduce  the  com- 
puter time  that  would  be  necessary  to  find  the  required  size  of 
the  slip  line  field  for  a preset  number  of  j lines. 

Initial  application  of  the  same  technique  to  the  present  prob- 
I’fem  created  difficulties,  because  the  number  of  j lines  may 
cliahge  in  the  various  slip  line  fields  when  inertia  forces  are 
taken  into  account.  In  the  iteration  technique  where  the  slip 
line  field  geometry  is  updated  for  velocity  increments  a variable 
number  of  j lines  is  undesirable,  because  for  an  updated  slip 
line  field  geometry  containing  more  j lines  than  the  previous 
field,  accelerations  for  points  on  the  additional  j line  are  not 
available.  For  this  reason  the  technique  using  a variable  number 
of  j lines  was  discontinued  in  this  program  and  a constant 
30  X 10  grid  size  was  adapted.  It  was  found  that  although  the 
computer  time  needed  to  perform  the  necessary  iterations  to  find 
the  required  size  of  the  slip  line  field  increased  considerably 
for  the  first  slip  line  field,  in  the  subsequent  iterative  proce- 
dure for  the  updating  of  the  slip  line  field  geometries  the 
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changes  in  the  size  of  these  fields  were  small  and  the  total  com- 
puter time  was  comparable  to  that  needed  with  variable  j lines. 

The  boundary  conditions  for  the  stress  field  also  presented 
problems  when  inertia  forces  were  included.  The  soil  surface  out- 
side the  tire-soil  contact  area  was  initially  assumed  as  stress 
free  and  the  direction  of  the  major  principal  stress  as  horizontal 
even  if  soil  inertia  forces  were  included.  Examination  of  the 
geometries  of  slip  line  fields  generated  with  this  assumption 
clearly  showed  an  abrupt  change  in  the  directions  of  slip  lines 
at  the  first  row  of  nodal  points  beneath  the  free  surface.  In 
further  analyses  the  direction  of  the  major  principal  stress  at 
the  surface  was  assumed  to  correspond  to  a hypothetical  loading  in 
the  direction  of  the  resultant  inertial  force.  Slip  line  fields 
generated  with  this  assumption  were  smooth  and  this  assumption  was 
adapted  in  the  further  development. 

The  singular  point  (Point  "A”  in  Fig.  2)  presents  special 
problems  when  inertial  accelerations  are  included  in  the  computa- 
tions. Stresses  at  this  point  are  computed  by  the  equation 

2 tan  cp(0-0  ) 

a = a e ° (16) 

that  is  the  solution  of  differential  equations  (1)  when  both  dx 
and  dz  vanish.  It  is  seen  that  the  inclusion  of  inertial  accel- 
erations does  not  affect  the  validity  of  Eq.  (16)  since  all  in- 
ertial terms  are  multiplied  by  the  vanishing  dx  or  dz  differ- 
ential. The  singular  point  can  be  regarded  as  a degenerate  j 
line  (with  zero  length)  along  which  9 changes  from  its  value  9^ 
at  the  adjacent  free  surface  to  that  specified  at  the  adjacent 

interface.  However,  for  the  determination  of  9 the  inertial 

o 

accelerations  would  have  to  be  known  at  the  singular  point  con- 
sidered as  part  of  the  free  surface.  For  the  computation  of  these 
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acceleration  velocities  at  this  point,  themselves  dependent  on  0 , 

o 

would  have  to  be  computed.  Since  the  neighborhood  of  the  singular 
point  and  the  singular  point  itself  is  a critical  part  of  the  solu- 
tion for  the  whole  slip  line  field,  it  was  necessary  to  develop 
and  adopt  an  updating  procedure  for  the  computation  of  6^  and 
velocities  for  this  point  independently  of  the  updating  of  the 
rest  of  the  field.  The  problem  of  updating  the  solution  of  the 
differential  equations  in  the  neighborhood  of  the  singular  point 
was  compounded  by  the  difficulties  of  applying  the  numerical  techr 
niques  to  the  solution.  The  assumption  of  a linear  variation  of 
velocities  from  an  i,j  point  to  an  i,j+l  and  i+l,j  point 
(applied  for  acceleration  computations  at  an  inner  i,j  point  of 
the  velocity  field)  is  clearly  inapplicable  to  the  singular  point 
that  is  the  common  location  of  nodal  points  with  various  i desig- 
nations. The  assumption  of  linearity  would  lead  in  this  case  to 
either  infinite,  or  undetermined  accelerations.  Thus,  it  was 
necessary  to  compute  the  theoretical  accelerations  at  the  singular 
point  from  another  combination  of  points  along  the  next  j line. 

An  extensive  analysis  of  the  acceleration  computation  was  made 
using  various  grid  sizes  and  combinations  of  two  or  more  points 
along  the  j+1  line.  From  these  analyses  it  was  concluded  that 
while  different  numerical  techniques  affected  the  value  of  accel- 
eration computed  for  this  point  the  problem  is  a physical  one: 
with  increasing  translational  velocity  the  location  where  acceler- 
ations become  critical  first  is  that  of  the  singular  point.  This 
could  be  expected  since  the  singular  point  is  also  the  entry  point 
where  the  tire  "hits'*  the  ground.  Accelerations  at  the  singular 
point  are  critical  when  the  inclination  of  the  major  principal 
stress  associated  with  these  accelerations  equals  that  at  which 
the  free  surface  becomes  a slip  line.  There  is  no  single  valued 
solution  for  the  slip  line  field  when  the  accelerations  at  this 
point  exceed  this  limit. 
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The  geometries  of  updated  slip  line  fields  generated  with  the 
iterative  procedure  for  the  computation  of  inertial  accelerations 
were  transmitted  to  a plotting  program  subroutine  that  allowed  the 
visual  examination  of  the  updated  geometries.  This  examination 
revealed  that  in  many  cases  when  the  numerical  end  results  appeared 
to  be  acceptable  the  numerical  solution  of  the  differential  equa- 
tions (1)  was  multivalued,  as  evidenced  by  the  overlapping  of  the 
slip  lines.  A slip  line  field  with  an  appreciable  overlap  is 
shown  in  Fig.  7.  In  many  cases  the  overlap  was  hardly  noticeable 


Fig.  7.  Overlap  of  Slip  Lines  in  the  Case  of  a Multivalued  Solution 
of  the  Governing  Differential  Equations 

upon  visual  examination  of  the  slip  line  field,  therefore,  a pro- 
vision was  made  in  the  program  to  indicate  by  a printout  if,  within 
the  accuracy  of  calculations,  overlap  occurs.  Overlap  indicates 
that  the  mathematical  solution  of  the  differential  equations  (1) 
at  a particular  location  is  multivalued,  i.e.,  more  than  one  a 
and  6 value  belongs  to  the  same  x,z  point  in  the  slip  line 
field. 

The  importance  of  an  overlap  indicating  multivaluedness  of 
the  mathematical  solution  is  that  the  solution  no  longer  represents 
the  physical  behavior  of  soil  under  the  applied  loads,  since  it  is 
a physical  impossibility  to  have  two  different  stress  states  in  the 
soil  at  the  same  place  and  at  the  same  time.  Thus,  overlap  is  a 
limitation  to  the  application  of  plasticity  theory  with  inertial 
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forces  to  the  tire-soil  interaction  problem.  Since  with  the  method 
of  velocity  fields,  overlapping  starts  to  occur  at  relatively  low 
velocities  (in  the  range  of  5 feet/second) ; this  limitation  is  a 
serious  one.  Therefore,  an  extensive  investigation  was  made  to 
determine  the  causes  of  the  overlap  and  especially,  whether  the 
overlap  is  the  result  of  the  application  of  numerical  techniques 
or  is  inherent  in  the  problem.  As  a result  of  these  investigations 
the  following  conclusions  were  drawn: 

1.  Overlap  occurs  if  there  is  an  abrupt  change  in 
the  magnitude  of  inertia  forces  between  adjacent 
nodal  points.  The  geometry  of  the  slip  line 
fields  is  such  that  while  there  are  no  abrupt 
changes  in  the  velocities  along  the  slip  lines, 
there  are  abrupt  changes  in  the  derivatives  of 
the  velocities  at  locations  where  the  various 
zones  (active,  passive,  and  radial)  are  joined 
together.  This  is  evident  in  the  case  of 
Prandtl  fields  for  weightless  soil  where  straight 
slip  lines  are  joined  with  logarithmic  spirals 

at  the  zone  boundaries.  Numerical  computation 
techniques  can  aggravate  this  problem  if  at  the 
zone  boundaries  a different  set  of  nodal  points 
is  used  for  the  acceleration  computation  than 
inside  the  field.  Numerical  computation  tech- 
niques can  alleviate  the  problem  by  using  schemes 
that  tend  to  smooth  these  abrupt  changes. 

2.  Overlap  usually  occurs  first  in  the  neighborhood 
of  the  singular  point.  Bow  waves  change  the 
geometry  of  free  surface  in  a way  that  the  ef- 
fect of  inertia  forces  in  the  neighborhood  of 
the  singular  point  is  counterbalanced. 
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3.  The  velocities  determined  by  the  theory  of 

velocity  fields  are  based  on  a two  dimensional 
model  of  the  tire-soil  interaction.  Implicit 
in  the  model  is  that  soil  particles  may  move 
only  in  the  plane  of  travel.  In  actuality  the 
movement  of  soil  particles  is  not  restricted 
to  the  plane  of  travel  and,  as  a consequence, 
the  actual  velocities  and  accelerations  are 
less  than  predicted  by  the  two  dimensional 
theory. 

Results  of  Sample  Computations 

For  the  evaluation  of  the  effect  of  inertia  forces  on  tire 
performance  sample  computations  were  made  within  the  velocity 
limits  imposed  by  the  criterion  that  the  solution  should  be  single 
valued.  The  results  of  these  sample  computations  were  surprising 
in  that  the  effect  of  inertia  forces  was  minimal  even  though  the 
computed  inertial  accelerations  often  exceeded  g by  a factor  of 
10  or  more.  Figure  8 shows  a comparison  of  normal  stresses  com- 
puted with  and  without  inertia  forces  generated  at  a translational 
velocity  of  5 feet /second.  There  is  a slight,  almost  impercepti- 
ble rise  in  the  normal  stresses  in  the  front  and  a very  slight 
decrease  of  them  in  the  rear  of  the  tire  when  inertia  forces  are 
accounted  for.  Table  1 shows  the  pertinent  input  and  performance 
data. 

Figure  9 shows  the  geometries  of  front  slip  line  fields  de- 
termined for  0,  3,  and  6 feet/second  translational  velocities, 
respectively.  Although  changes  in  the  geometry  with  the  velocity 
appear  to  be  minor,  closer  observations  of  the  "i"  lines  shows 
their  approaching  each  other  in  the  vicinity  of  the  boundary  be- 
tween the  passive  and  radial  zone,  an  indication  of  an  impending 
overlap  situation. 
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Pig.  8.  Normal  Stress  Distribution  V/ith/V/ithout  Soil  Inertia  Forces 
for  Tire  Loading  and  Soil  Conditions  Shown  in  Table  1. 

Summary  Discussion  of  the  Method  of 
Velocity  Fields  and  Conclusions 

A computer  program  has  been  developed  that  takes  soil  inertia 
forces  in  tire-soil  interaction  into  account  by  computing  the 
velocity  fields  associated  with  the  slip  line  fields  and  the  ac- 
celerations therefrom.  It  was  found  that  the  method  can  be  used 
up  to  about  5 feet/second  velocity  with  9.00-14  tires  (for 
larger  tires  probably  up  to  proportionately  higher  velocities) 
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Fig.  9.  Geometry  of  Front  Slip  Line  Fields  Computed  For  v- 
(No  Soil  Inertia  Forces)  and  = 3 and  6 ft/sec  ^ 
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TABLE  1 TIRE  PERFORMANCE  WITH  AND  WITHOUT 
SOIL  INERTIA  FORCES 

(Method  of  Velocity  Fields) 


Tire  Size: 

Tire  Radius  (nominal) : 
Tire  Width: 

Inflation  Pressure  (Pj|^): 
Slip: 

Angle  of  Interface 
Friction: 


9.00-14 
1.18  ft 
0.70  ft 
12.5  psi 
15% 

15. 7^^ 


Front  Field 


Rear  Field 


Cohesion 
Friction  Angle 
Unit  Weight 

Translational  Velocity 
Load 

Drawbar  Pull 

Torque 

Sinkage 

Pull  Coefficient 


15  Ibs/sq  ft 
24° 

100  Ibs/cu  ft 
= 0 ft/sec 
408  lbs 
48  lbs 
129  ft-lb 
2.03  in. 

0.119 


30  Ibs/sq  ft 
28° 

105  Ibs/cu  ft 
= 5 ft/sec 
410  lbs 
53  lbs 
135  ft-lb 
1.81  in. 

0.129 


when  the  solution  of  the  differential  equations  becomes  multivalued. 
Because  of  the  two  dimensionality  of  the  velocity  fields,  veloci- 
ties and  associated  inertial  accelerations  are  overestimated  by 
this  method.  It  was  found  that  abrupt  changes  in  inertia  forces 
occur  inherently  in  this  method  because  derivatives  of  slip  line 
directions  that  control  the  velocity  field  are  discontinuous  at 
the  joint  boundaries  of  the  active,  passive,  and  radial  zones. 
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There  is  no  experimental  information  on  the  physical  phenome- 
non that  occurs  when  the  plasticity  solution  becomes  multivalued. 
Speculatively,  one  may  assume  that  failure  in  the  classical  Coulomb 
sense  cannot  occur  and  a so-called  "rigid"  soil  zone  exists  within 
the  bounds  of  overlap.  It  is  also  likely  that  the  soil  adjusts 
itself  and  seeks  to  fail  along  slip  lines  with  continuous  curva- 
ture. Upper  bound  solutions  based  on  the  characteristics  of  the 
differential  equations  for  velocities  would  have  to  be  considered 
rather  than  solutions  based  on  stress  characteristics.  A possible 
solution  would  be  a slip  line  field  with  only  one  zone  and  con- 
tinuous curvature  of  the  slip  lines. 

Inertia  forces  appear  to  have  little  effect  on  tire  perfor- 
mance as  long  as  the  plasticity  solution  is  single  valued.  A 
slight  rise  in  the  normal  stresses  in  the  front  field  increases 
tire  deflection  thus  resulting  in  a slightly  better  tractive  per- 
formance. It  should  be  noted  that  even  though  the  calculations 
were  made  for  a low  translational  velocity  the  calculated  inertia 
forces  inside  the  field  were  in  the  range  of  several,  often  more 
than  10  g-s. 

An  incidental  but  important  result  with  this  method  is  the 
introduction  of  a new  concept,  that  of  "contact  slip"  in  the  tire 
soil  interaction  problem.  For  determination  of  the  contact  slip 
in  any  interaction  problem,  only  a small  portion  of  the  computer 
program,  dealing  with  the  kinematic  boundary  conditions,  is  needed. 
The  new  concept  can  be  the  basis  of  further  theoretical  research 
in  the  area  of  slip,  mobilized  shear,  and  traction  development. 


V.  DETERMINATION  OF  SOIL  INERTIA  FORCES 
BY  THE  PARTICLE  PATH  METHOD 


Introduction 

Soil  inertia  forces  are  generated  by  the  acceleration  and 
deceleration  of  soil  particles  during  the  passage  of  a tire.  Ihe 
motion  of  soil  particles  follows  a particle  path  that  defines  the 
geometry  of  particle  motion  in  a coordinate  system  fixed  to  the 
ground.  Velocities  and  accelerations  are  time  derivatives  of  the 
coordinate  vectors  of  particle  path. 

Experimental  Information  on  Particle  Path  Geometry 

Measurements  of  particle  paths  beneath  rigid  wheels  were  per- 
formed by  several  researchers  under  laboratory  conditions  (Refs.  20 
through  22)  using  the  flash  X-ray  technique.  These  experiments 
show  the  same  general  pattern  and  dependence  of  the  particle  path 
on  slip,  as  shown,  for  example,  in  Fig.  10.  Although  not  pointed 
out  in  the  referenced  publications  that  describe  these  particle 
path  measurements,  theoretical  considerations  indicate  that  for  a 
constant  velocity  the  particle  path  geometry  must  not  change  with 
the  X coordinate,  the  direction  of  travel.  Any  variation  in 
particle  path  geometry  in  the  x direction  obtained  in  experiments 
must  be  the  result  of  experimental  error,  inhomogeneities  in  the 
soil  bed  preparation,  etc.  The  invariability  of  the  particle  path 
geometry  with  respect  to  the  direction  of  travel  also  follows  from 
the  "steady  state"  concept  of  tire-soil  interaction. 

Hie  experiments  show,  on  the  other  hand,  that  the  size  of 
particle  paths  diminishes  with  the  depth  beneath  the  wheel  and 
there  is  a certain  limit  depth  below  which  the  soil  is  unaffected  ^ 

by  the  passage  of  the  wheel.  I 
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Fig.  10.  Particle  Motion  as  Influenced  by  Slip  Rate 
(From  Ref.  21) 

There  are  some  conditions  that  define  the  end  points  of  the 
path  of  particles  that  were  originally  at  the  surface.  After  the 
passage  of  the  tire  these  particles  wind  up  in  the  rut  at  some 
depth,  that  evidently  equals  the  vertical  distance  between  the  end 
points  of  the  path.  Experiments  also  show  that  the  horizontal  posi- 
tion of  a particle  at  the  surface  is  changed  by  the  action  of  the 
wheel  as  it  passes  over  the  particle.  In  the  case  of  a slipping 
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wheel,  the  particle  at  the  surface  is  pushed  backward  from  its 
original  position;  in  the  case  of  a skidding  wheel,  or  negative 
slip,  the  particle  at  the  surface  is  pushed  in  the  forward  direc- 


tion. Figure  11  shows  the  displacement  of  dyed  sand  particles 


Fig.  11.  Displacement  of  Sand  Particles  Under  a Slipping  Wheel 
(a,  s = 55^)  and  a Skidding  Wheel  (b,  s = 


beneath  a skidding  and  a slipping  wheel  determined  in  the  mobility 
bin  at  Grumman.  Thus,  a particle  at  tiie  surface  is  generally  dis- 
placed both  in  the  horizontal  and  vertical  direction;  the  displace- 
ment equals  the  difference  between  the  origin  and  the  terminus  of 
the  path  of  a particle  at  the  surface. 

Idealization  of  Particle  Path  Geometry 

For  the  intended  use  of  particle  path  geometries  it  is  desir- 
able to  describe  these  geometries  by  analytical  functions.  Experi- 
mentally determined  particle  paths  have  some  resemblance  to 
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cardioids  and  in  some  publications  they  are  referred  to  as  such. 
The  parametric  equation  of  a cardioid  is  (Fig.  12) 


X = 2a  sin  o)  - a sin  2co 

z *=  2a  cos  0)  ~ a cos  2cd 
z 


(17) 


Fig.  12  Cardioid  Geometry 

Obviously,  a cardioid  is  not  a suitable  curve  to  describe  the 
parcicle  path  when  constraints  are  set  for  the  displacement  of  the 
end  point  in  both  the  horizontal  and  vertical  direction. 

Another  curve  that  has  a shape  similar  to  the  observed  parti- 
cle paths  is  the  nephroid.  The  parametric  equation  is  (Fig.  13) 

X = 3a  sin  o)  - a sin  3o> 

, (18) 

z = 3a  cos  o>  - a cos  3cu 
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Equation  (19)  refers  to  a coordinate  system  with  its  origin 
coinciding  with  the  initial  at  rest  position  of  the  particle.  The 
expansion  of  the  expression  in  parentheses  for  z by  the  term 
(p2  ■ 1)  serves  this  purpose. 

Equation  (19)  contains  parametric  exnressions  of  the  particle 
path  that  represents  a family  of  curves  of  various  shapes.  Fig- 
ure 14  shows  variations  of  the  geometry  as  defined  by  Eq.  (19)  for 
various  values  of  the  parameters  and  P2*  Equation  (19)  also 

represents  the  geometry  of  the  particle  path  in  parametric  form. 
Time  derivatives  of  the  coordinates  x and  z yield  the  veloci- 
ties and  accelerations  of  the  particle  at  the  time  the  particle  is 
at  point  x,z.  To  calculate  these  time  derivatives  the  parameter 
d has  to  be  related  to  time,  t,  preferably  by  a differentiable 
function.  If 

0)  = f(t)  (20; 

then  the  time  derivatives  of  the  coordinates  of  the  particle  path 
are 

• • 

X = ap^o)(cos  o>  - cos  p^^co) 

• • 

z = ap2<J^(sin  od  - sin  P2CD) 

0 (21 

••  ••  ♦ y ' ' 

X = ap^co(cos  0)  - cos  p^oo)  + ap^(co)  (-sin  a>  + p^  sin  p^^o)) 

••  ••  *2 

z = ap2^i>(sin  co  - sin  P2«^)  + ap2(a>)  (cos  cos  P2O0) 

To  compute  the  accelerations  from  Eq.  (21)  the  relationship 
between  the  parameter  o)  and  time,  t,  has  to  be  defined.  To 
this  end  it  is  useful  to  consider  the  changes  in  the  direction  of 
the  particle  movement  that  can  be  associated  with  the  relative  po- 
sition of  the  tire.  The  direction  of  the  particle  at  any  point 


along  the  particle  path  is  given  by  the  tangent  of  the  path  at 
that  point.  From  Eq.  (19)  the  tangent  of  the  particle  path  can  be 
calculated  as 


dz 

dz  _ ^ 
dx  ^ 

do) 


P2  sin  CO  - sin 
cos  CO  - cos  p^co 


(22) 


The  points  along  the  particle  path  that  may  be  associated  with 
various  positions  of  the  tire  are  shown  in  Fig.  15.  The  value  of 
the  parameter  at  these  points  may  be  determined  from  the  condition 
that  the  tangent  of  the  particle  path  is  either  horizontal  or  ver- 
tical. The  following  tabulation  shows  the  co  values  at  these 
points. 


TABLE  2 VALUES  OF  THE  PARAMETER  co 
AT  VARIOUS  POINTS  OF  PARTICLE  PATH 


Point  Tangent  co 


A 

B 

C 

D 

E 


Vertical 

Horizontal 

Vertical 

Horizontal 

Vertical 


0 

7t/(1  + p^) 
27t/  (1  + p^) 
3tt/(1  + P2) 
2tt/(p2  - 1) 


The  position  of  the  tire  that  may  be  associated  with  the 
above  points  and  the  time  elapsed  while  the  tire  passes  from  one 
point  to  another  is  shown  in  the  lower  part  of  Fig.  15.  Obviously, 
there  is  no  rigorous  time  relation  between  these  particle  locations 
and  tire  travel  positions  and  these  considerations  serve  primarily 
as  a guide  for  the  mathematical  formulation  of  the  cd  = f(t)  re- 
lationship. 
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As  mentioned  before,  for  a particle  at  the  surface  the  end 
points  also  have  to  meet  certain  conditions.  First,  the  vertical 
distance  between  the  end  points  must  equal  the  rut  depth,  i.e.. 


z 


E 


(23) 


The  constant  "a"  in  Eq.  (19)  is  selected  so  that  this  condition 
be  met. 


Second,  the  assumption  is  made  that  the  horizontal  distance 
between  the  end  points  equals  the  contact  slip  times  the  rut  depth, 
or 


If  in  Eq.  (19)  the  x,z  values  are  substituted  with  the 
values  given  by  Eqs.  (23)  and  (24)  and  with  the  appropriate  value 
of  the  o)  parameter  given  in  Table  2 then  a relationship  between 
and  s'  can  be  obtained.  This  relationship  is  shown  in  Fig- 
ure 16  for  p2  values  varying  from  0.75  p^  to  0.95  p^.  It  is 
seen  that  within  these  limits  the  value  of  Pj^  is  affected  but 
little  by  the  value  of  P2»  and  Fig.  16  could  be  used  for  the  es- 
timation of  p^  for  any  value  of  the  contact  slip  s'.  On  the 
basis  of  the  available  experimental  information  on  particle  paths, 
it  appears  that  a P2  value  of  about  0.9  Pj^  is  a reasonable 
assumption.  Should  further  information  on  particle  paths  become 
available,  any  other  value  of  P2  that  allows  a good  simulation 
of  particle  path  may  be  selected.  In  this  case  the  value  of  p^ 
may  be  computed  from  the  following  pol3nnomial  determined  as  best 
fitting  the  curves  shown  in  Fig.  16. 
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where 


3 2 


+ a 


33  = - 717.661  + 2504.05  - 3253.0 

32  = - 10.986  + 24.6844  ^ - VJ .ei9h  i + 4.71601 

3j^  = .273714  + 1.05039  i - 2.0157 

3q  = - 1.40114  + 2.37218  ^ + 2.01924 


The  parameters  Pi>P2  ^ completely  define  the  particle 

path  and  the  accelerations  if  to  = f(t)  is  known.  Various  func- 
tional forms  for  this  relationship  have  been  assumed  and  investi- 
gated for  suitability  to  fit  the  time  constraints  established  pre- 
viously for  the  various  points  along  the  path.  Polynomials  gen- 
erally did  not  meet  the  requirement  that  cd  should  increase 
monotonically  with  t.  An  analytical  study  of  various  logarithmic 
and  exponential  expressions  showed  that  if  the  requirement  for  a 
monotonic  increase  of  00  with  t was  satisfied  then,  in  some 
instances,  it  was  not  possible  to  meet  the  time  constraints. 
Finally,  the  following  expression  was  chosen 

00  = Cj^  arctan  (c2(t  + + co^  (26) 

The  constants  c^,C2  in  Eq.  (26)  are  determined  from  the 
ratios  of  the  elapsed  time  between  particle  path  locations  A,B,D, 
and  E.  The  value  of  00^  in  Eq.  (25)  may  be  tentatively  chosen  as 
0.5  X 0)  . 

Cl 
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Use  of  the  Idealized  Particle  Path  Geometry 
in  the  Tire-Soil  Model 


The  anal3rtical  expression  for  the  particle  path  geometry  is 
convenient  for  the  determination  of  accelerations.  To  use  this 
analytical  expression  for  the  calculation  of  inertia  forces  in  the 
tire-soil  model  it  is  necessary  to  determine  the  constants  in  that 
expression.  To  this  end  it  is  assumed  that  the  constants  in  the 
particle  path  equation  may  be  determined  from  the  * 0 case. 

An  alternative  to  this  assumption  is  that  the  particle  path  re- 
mains unchanged  for  small  increments  of  velocity  and  the  particle 
path  geometry  for  a given  translational  velocity  would  be  deter- 
mined by  an  updating  procedure.  Analysis  of  a few  cases  indicated 
that,  in  view  of  the  approximate  nature  of  the  particle  path  con- 
cept, an  elaborate  updating  procedure  would  be  justified  only  if 
more  experimental  information  on  the  geometry  of  particle  paths 
were  available. 

For  the  assumption  chat  the  particle  path  geometry  may  be  de- 
termined from  the  v^  = 0 case,  computations  of  accelerations 
become  a simple  matter.  They  require  the  addition  of  a short  sub- 
routine named  "ACCE"  to  the  program  that  was  developed  for  the 
computation  of  tire  performance  for  given  tire  and  soil  input  data. 
(Main  program  "KTIRE"  and  subroutine  "SLFI,"  developed  for  tire- 
soil  interaction  simulation,  are  described  in  TACOM  Tech.  Report 
No.  11900  (LL147)  (Ref.  4).  The  subroutine  "ACCE"  is  called  from 
subroutine  "SLFI"  after  the  first  approximation  of  the  coordinates 
x,z  is  computed  for  each  nodal  point  of  the  slip  line  field.  The 
accelerations  computed  by  the  subroutine  "ACCE"  are  then  used  in 
the  difference  equations  for  the  computation  of  a and  d and 
for  the  computation  of  the  improved  values  of  the  x,z  coordinates 
of  the  nodal  point. 
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A detailed  flow  chart  of  the  expanded  program  for  the  computa- 
tion of  tire  performance  with  consideration  of  soil  inertia  force 
is  given  in  Appendix  B. 

Problems  Encountered  with  the  Development  of 
the  Computer  Program 

Ihese  were  generally  of  the  same  nature  as  with  the  program 
based  on  velocity  fields.  However,  the  analytical  simulation  of 
the  particle  path  geometry  resulted  in  continuous  time  derivatives 
alleviating  the  overlap  problems  at  zone  transitions  encountered 
with  the  other  program.  The  critical  accelerations  occurring  at 
the  singular  point,  that  are  inherent  in  the  problem,  limit  the 
applicability  of  the  method  to  about  10  to  20  feet/second  trans- 
lational velocity,  depending  on  the  tire  dimensions. 

In  cases  where  the  soil  and  loading  conditions  are  such  that 
there  is  no  slip  line  field  in  the  front,  association  of  the  vari- 
ous points  of  the  particle  path  with  relative  positions  of  the 
tire  is  no  longer  applicable.  Experimental  information  on  parti- 
cle paths  generated  under  such  conditions  is  needed  for  meaningful 
simulation  of  the  geometry  of  particle  paths  and  the  motion  of 
particles  in  the  time  frame  of  tire  passage. 

Results  of  Sample  Computations 

For  evaluation  of  the  effect  of  soil  inertia  forces  on  tire 
performance  by  this  method,  computations  were  performed  for  a few 
typical  cases.  Results  of  these  computations  are  shown  in  Tables  3 
and  4. 

In  the  computations  for  the  performance  calculations  shown  in 
Tables  3 and  4 no  overlap  was  experienced  for  the  translational 
velocities  indicated,  i.e.,  all  slip  line  field  solutions  were 
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TABLE  3 TIRE  PERFORMANCE  AT  VARIOUS  VELOCITIES 


Tire  Size: 

9.00-14 

Tire  Radius: 

R « 1.18 

ft 

Tire  Width: 

B * 0.74 

ft 

Load: 

L = 614 

lbs 

Inflation  Pressure: 

P£  “ 10. 

9 psi 

Limit  Pressure: 

Pi  - 11-' 

0 psi 

Deflection  Coefficient: 

€ » 0.9 

Slip: 

s = 10% 

Cone  Index  Gradient: 

CGR  = 15 

pci 

Interface  Friction  Angle: 

6 = 22. 

8° 

Translational 

Pull 

Sinkage 

Contact 

Entry 

Rear 

Velocity  ft/sec  Coefficient  in. 

Slip 

Angle  (°) 

Angle  ( 

0 

0.3268 

0.724 

0.0738 

26.06 

12 

5 

0.3334 

0.599 

0.0614 

23.90 

10 

8 

0.3340 

0.593 

0.0600 

23.84 

10 

10 

0.3346 

0.588 

0.0618 

23.77 

10 

12 

0.3352 

0.583 

NA 

23.71 

10 

single  valued. 

The  highest  velocities  shown 

in  these  tables  are 

approximately  the  limits  above  which  the  solution  becomes  multi- 
valued . 


It  is  seen  that  the  increase  of  tire  performance  due  to  soil 
inertia  forces  are  minimal  in  these  cases  and,  for  all  practical 
purposes,  negligible. 


49 


TABLE  4 TIRE  PERFORMANCE  AT  VARIOUS  VELOCITIES 


Tire  Size: 

4.00-7 

Tire  Radius: 

R = 0.59 

ft 

Tire  Width: 

B » 0.40 

ft 

Load: 

L = 225  lbs 

Inflation  Pressure: 

p.  = 10.3 

psi 

Limit  Pressure: 

= 10.6 

psi 

Deflection  Coefficient: 

€ = 0.85 

Slip: 

s = 30% 

Cone  Index  Gradient: 

CCR=  19.6 

pci 

Interface  Friction  Angle: 

6 = 27.8 

0 

Translational  Pull 

Sinkage 

Contact 

Entry 

Rear 

Velocity  ft/sec  Coefficient  in. 

Slip 

Angle  C) 

Angle 

0 0.3217 

1.11 

0.089 

42.3 

22 

5 0.3270 

1.08 

0.093 

42.0 

22 

8 0.3291 

1.07 

0.131 

41.9 

22 

10  0.3508 

0.93 

0.135 

39.3 

20 

Summary  Discussion  of  the  Particle  Path  Method 
and  Conclusions 

Uie  consideration  of  soil  inertia  forces  by  the  particle  path 
method  is  simpler  than  with  the  method  of  velocity  fields.  The 
analytical  simulation  of  particle  path  geometries  results  in  smooth 
changes  of  velocities  and  accelerations  from  one  nodal  point  to 
another  in  the  slip  line  field.  As  a result,  the  limitation  im- 
posed by  the  raultivaluedness  of  the  solution  of  the  differential 
equations  defining  the  slip  line  fields  is  less  severe  than  the 
velocity  field  method,  and  analysis  of  tire  performance  at  higher 
velocities  becomes  possible. 
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The  drawback  of  the  particle  path  method  is  that  its  enq)iri- 
cal  basis  includes  only  a few  tests.  Therefore,  it  was  unavoidable 
to  make  somewhat  arbitrary  assumptions  in  the  analytical  simulation 
of  the  particle  path  geometry  to  attain  a usable  method.  However, 
there  is  enough  freedom  in  the  anal3rtical  simulation  of  particle 
paths  to  allow  for  adjustments  should  more  experimental  information 
become  available. 

The  most  interesting  and  important  result  of  the  analyses  per- 
formed by  the  method  of  particle  paths  is  that,  in  agreement  with 
the  results  obtained  with  the  issthod  of  velocity  fields,  the  effect 
of  soil  inertia  forces  on  tire  performance  is  minimal  for  the 
translational  velocities  investigated.  An  explanation  for  this 
finding,  that  is  somewhat  contrary  to  expectations,  is  that  in  the 
cases  analyzed  tire  performance  was  governed  by  failure  conditions 
in  the  soil.  Even  though  soil  inertia  forces  affected  the  geometry 
of  slip  line  fields,  the  stresses  in  the  soil  were  controlled  by 
its  strength;  therefore,  soil  inertia  forces  had  little  effect  on 
the  interface  stresses  and,  consequently,  on  tire  performance.  A 
logical  amplification  of  this  explanation  is  that  soil  inertia 
forces  are  very  likely  to  have  appreciable  effect  on  tire  perfor- 
mance if  they  affect  the  mode  of  failure  in  the  soil.  Multivalued- 
ness of  the  solution  of  Eq.  (1)  is  an  indication  that  the  soil  may 
not  fail  according  to  the  concepts  of  soil  plasticity.  At  this 
time  neither  theory  nor  experiments  are  available  that  would  apply 
to  the  problem. 

Another  result  of  the  sample  calculations  worth  noting  is 
that  the  computed  average  value  of  contact  slip  changes  with  veloc- 
ity. The  improvements  in  tire  performance  with  speed  reported  in 
Ref.  1 refer  to  a constant  slip  of  20  percent.  The  question  may 
be  raised  whether  the  improvement  in  tire  performance  with  speed 
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would  manifest  itself  over  the  entire  slip  range,  including  the 
maximum  pull  performance,  or  whether  the  observed  improvements 
signify  only  an  accelerated  development  of  the  interface  shear 
stresses  with  slip  as  the  travel  velocity  increased.  Tire  perfor- 
mance experiments  at  various  velocities  over  the  entire  slip  range 
are  needed  to  clarify  this  problem. 
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VI.  EFFECT  OF  LOADING  RATE  ON  SOIL  STRENGTH 
AND  TIRE-SOIL  INTERACTION 


Introduction 

In  tire-soil  interaction  the  soil  is  loaded  at  a rate  that  is 
directly  proportional  to  the  translational  velocity  of  the  tire. 

The  rate  of  loading  affects  the  strength  of  soil  in  many  ways.  In 
civil  soils  engineering  practice  consolidation  and  creep  theories 
deal  with  the  long  term  effects  of  loading  rates  on  soil  strength; 
the  l.oading  rates  encountered  in  tire-soil  interaction  problems 
fall  into  the  category  termed  as  "rapid"  or  "quick"  loading  and  the 
soil  strength  referred  to  as  "undrained"  indicating  that  during 
such  rapid  loading  there  is  no  drainage  of  the  soil.  For  the  pur- 
poses of  civil  soils  engineering,  laboratory  tests  aimed  at  deter- 
mining the  strength  of  soil  under  rapid  loading  conditions  are 
performed  so  as  to  reach  the  maximum  load  in  a matter  of  minutes. 

In  tire-soil  interaction  the  loading  rate  is  generally  faster 
than  that  considered  "rapid"  in  civil  soil  engineering.  Soil  be- 
havior at  loading  rates  comparable  to  those  encountered  in  tire- 
soil  interaction  has  not  been  investigated  until  recently,  since 
laboratory  investigations  at  these  loading  rates  require  special 
equipment,  referred  to  as  "dynamic"  triaxial  testing  apparatus. 

In  dynamic  tests,  loads  are  applied  in  fractions  of  seconds.  Dy- 
namic tests  are  generally  strain  controlled  and  investigations 
relate  strength  properties  to  the  strain  rate. 

Indirect  information  on  the  effect  of  loading  rate  on  the 
strength  properties  of  soils  is  also  available  from  various  sources 
where  the  resistance  of  soil  to  the  penetration  of  plates,  blades, 
or  cones  is  measured  at  various  rates  of  penetration. 
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Physical  Causes  of  Strain  Rate  Dependent 
Strength  Properties  of  Soils 

The  physical  causes  of  the  strain  rate  dependent  behavior  of 
soils  are  manifold  and  some  of  them  are  not  clearly  understood. 

The  rate  of  strain  affects  soil  behavior  on  the  microlevel  through 
its  effect  on  the  physicochemical  and  electrochemical  bonding 
forces  between  the  soil  particles.  On  the  macrolevel,  moisture 
redistribution  during  shearing  and  plastic  flow  of  contacting  par- 
ticle asperities  are  the  major  time-dependent  phenomena  that  affect 
soil  strength  and  tire-soil  interaction.  Since  these  phenomena 
affect  the  results  of  the  various  types  of  laboratory  and  field 
strength  tests  and  tire-soil  interaction  in  a different  way,  it  is 
important  to  consider  the  mechanism  underlying  these  phenomena 
when  strain  rate  effects  are  analyzed.  In  the  following  discussion 
many  aspects  of  these  very  complex  phenomena  are  omitted  so  that 
the  essential  mechanism  may  be  clearly  presented.  For  the  purposes 
of  this  discussion  it  is  convenient  to  group  soils  in  three  cate- 
gories: saturated  cohesive  type  soils,  cohesionless  soils,  and 

partially  saturated  frictional -cohesive  soils. 

In  saturated  cohesive  type  soils,  such  as  the  Buckshot  clay, 
moisture  redistribution  during  loading  is  the  most  important  single 
factor  responsible  for  the  increase  of  soil  strength  with  the  rate 
of  strain.  The  phenomenon  of  moisture  redistribution  can  be  best 
explained  by  considering  a triaxial  test  (Fig.  17) . Loading  of 
the  sample  in  a triaxial  test  results  in  a stress  state  ultimately 
leading  to  the  development  of  a failure  surface  where  shear 
stresses  exceed  the  shear  strength  of  the  material.  In  the  process 
of  straining  the  sample,  high  shear  stresses  along  the  potential 
failure  surface  tend  to  dilate  the  soil.  This  tendency  to  dilate 
reduces  the  pore  water  pressures  along  the  potential  failure 
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Fig.  17.  Moisture  Migration  in  Triaxlal  Test 

surface  relative  to  the  rest  of  the  sample.  The  reduced  pore  water 
pressures  suck  moisture  toward  the  failure  surface,  therefore  in- 
creasing the  moisture  content  locally.  Thus,  a moisture  redis- 
tribution occurs  in  the  sample  during  the  test  and  the  actual 
strength  of  the  sample  corresponds  to  the  moisture  content  in  the 
immediate  vicinity  of  the  failure  surface.  Experiments  confirm 
that  moisture  redistribution  takes  place  in  dynamically  tested 
triaxial  samples  (Ref.  23) . 

The  rate  of  moisture  migration  depends  on  the  permeability  of 
the  clay.  The  total  amount  of  moisture  that  migrates  during  a test 
depends  on  the  rate  of  migration  and  the  time  available  for  migra- 
tion. Obviously,  this  latter  depends  on  the  strain  rate  and  hence 
dependence  of  the  strength  on  the  strain  rate.  It  should  be  noted 
that  this  moisture  migration  occurs  even  though  the  test  is  of  the 
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"undrained"  t3rpe,  i.e.,  no  water  is  going  in  or  out  to  the  sample 
during  the  test  and  the  total  water  content  of  the  sample  has  not 
changed. 

The  process  of  moisture  redistribution  is  extremely  complex 
and  its  analjrtical  treatment  would  require  the  knowledge  of  stress 
states  in  the  specimen  for  the  deformed  geometry  as  well  as  pore 
water  pressure  response  of  the  soil  to  these  stress  states.  From 
the  viewpoint  of  tire-soil  interaction,  a significant  feature  of 
this  complex  redistribution  process  is  that  the  rate  of  moisture 
migration  toward  a potential  failure  surface  depends  not  only  on 
the  properties  of  soil,  but  also  on  the  stress  states  in  the 
vicinity  of  that  surface  and  the  distance  from  the  impervious 
boundaries  (the  sealed  bottom  and  the  rubber  membrane  in  the  tri- 
axial  test)  of  the  specimen.  These  conditions  are  different  in 
tire-soil  interaction  from  those  in  a triaxial  or  a field  test, 
therefore,  strength-strain  rate  relations  found  by  such  tests  are 
not  immediately  applicable  to  tire-soil  interaction  problems. 

In  cohesionless  soils  the  experimental  evidence  concerning 
the  effect  of  strain  rate  on  their  strength  is  inconclusive. 

Early  experiments  by  Casagrande  and  Shannon  (Ref.  24)  showed  a 
small  increase  in  strength  with  an  increase  of  the  strain  rate, 
while  Whitman  and  Healy  found  the  shear  strength  of  sands  unaf- 
fected by  the  strain  rate  (Ref.  25).  Larew  and  Atakol  (Ref.  26) 
found  that  the  rate  of  shear  strain  affected  the  strength  of  co- 
hesionless soils.  These  experimental  results  indicate  that  the 
response  of  cohesionless  soils  to  dynamic  loading  is  not  uniform 
and  the  mineral  composition  of  the  grains  and  arrangement  of  the 
particles  also  play  a role.  The  strain  rate  may  affect  this  re- 
sponse by  delaying  the  force  transmission  from  particle  to  parti- 
cle due  to  the  inertia  of  the  individual  particles  and  the  plastic 
deformation  of  the  contacting  surfaces. 
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In  partially  saturated  frictional-cohesive  soils  the  effect 
of  strain  rate  on  strength  properties  is  diverse.  In  the  three- 
phase  solid-water-air  system  pressures  of  pore  water,  as  well  as 
of  pore  air,  may  be  generated  and  moisture  redistribution  during 
loading  may  be  prevented  by  pore  air  pressures.  Also,  if  there 
are  open  channels,  the  pore  air  pressure  may  be  relieved;  the  time 
rate  of  this  relief  would  depend  on  the  air  permeability  of  soil. 

In  summary,  the  effect  of  strain  rate  on  the  strength  properties 
of  partially  saturated  frictional -cohesive  soils  is  exceedingly 
complex  and  cannot  be  evaluated,  even  qualitatively,  on  the  basis 
of  theoretical  and  experimental  information  available  at  present. 

Correlation  of  Strain-Rate  Effects  in  Laboratory 
and  Field  Tests  for  Strength  Property  Determination 
and  in  Tire-Soil  Interaction 

Dynamic  triaxial  tests  performed  under  controlled  laboratory 
conditions  are  the  best  source  of  information  on  the  relationship 
between  strain  rate  and  strength  in  various  tjrpes  of  soils.  Even 
so,  the  measured  resistance  of  the  soil  sample  to  the  applied 
loading  includes  the  effects  of  soil  inertia.  Only  in  the  most 
recent  experiments  (Ref.  27)  has  an  attempt  been  made  to  separate 
inertia  effects  in  the  triaxial  test  from  the  shearing  resistance. 

In  the  dynamic  triaxial  tests  the  applied  strain  rate  is  con- 
stant and  controlled  during  the  test.  Thus,  experimental  relation- 
ships between  strain  rate  and  strength  can  be  clearly  formulated; 
these  relationships  may  be  approximated  by  a logarithmic  relation- 
ship between  strain  rate  and  strength.  Figure  18,  from  Ref.  28, 
shows  such  a relationship  for  normally  and  overconsolidated  clays. 

Indirect  experimental  information  on  the  effect  of  strain 
rate  on  the  strength  properties  of  soils  has  also  been  obtained 
from  dynamic  plate  bearing  and  cone  penetration  tests  (Refs.  29 
through  31) . In  these  experiments  the  total  force  required  to 
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Fig.  l8.  Soil  Strength  - Strain  Rate  Relationships  for  Normally 
and  Overconsolidated  Clays  (Prom  Ref.  28) 

force  the  plate  or  cone  into  the  soil  was  measured  at  various  rates 
of  penetration.  While  these  experiments  provide  useful  information 
on  the  behavior  of  the  soils  tested  under  these  conditions  there  is 
no  theoretical  basis  for  the  generalization  of  these  findings  to 
other  types  of  soils.  Also,  the  evaluation  of  the  strength  prop- 
erties of  soils  from  these  tests  is  difficult  for  two  reasons. 
First,  the  effect  of  soil  inertia  forces  would  have  to  be  accounted 
for  separately.  Second,  the  strain  and  the  strain  rate  in  the  soil 
volume  affected  varies,  and,  therefore,  any  observed  effect  of  the 
penetration  rate  on  the  load  is  the  cumulative  effect  of  various 
strain  rates  in  the  soil. 


To  compare,  at  least  approximately,  results  of  cone  penetra- 
tion tests  and  dynamic  triaxial  tests,  an  average  strain  rate  was 
calculated  for  cone  penetration  tests  performed  in  Buckshot  clay 
on  the  following  basis.  Slip  line  fields  for  cone  penetration  in 
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clay  show  that  a soil  volume  of  a radius  of  about  R = (5  to  6)  * r 

(r  « radius  of  the  base  of  the  cone)  is  affected  by  the  penetration 

of  the  cone.  If  the  initial  void  ratio  of  the  soil  is  e and 

o 

R = 5r  is  assumed,  then  the  change  of  void  ratio  due  to  the  volume 
displaced  by  the  cone  (assuming  no  surface  heave)  is 


e^  - e + 1 
_ o o 

24(1  + e^) 


(27) 


The  corresponding  strain  is 


e 


Ae 

1 + e 

o 


(28) 


The  average  strain  rate  depends  on  the  average  strain  and  the  rate 
of  penetration.  For  the  WES  cone  the  strain  rate  is 


€ = 


(29) 


where  h ® 3.77  cm,  the  height  of  the  cone.  For  the  standard 
rate  of  cone  penetration  (3.05  cm/sec) 


e = .8  X €(sec  ^)  (30) 

In  the  tire  tests  conducted  in  Buckshot  clay  the  cone  index 
varied  from  16  to  62.  According  to  experiments  reported  in 
Ref.  32  this  cone  index  variation  corresponds  to  a variation  in 
moisture  content  from  40  to  33  percent.  The  corresponding  ini- 
tial void  ratios  vary  from  1.37  to  0.97,  the  average  strains 
from  1.2  to  1.0  percent,  and  the  average  strain  rates  from 
0.8  to  1.0  percent  sec 

In  tire-soil  interaction  a crude  approximation  for  the  average 
strain  ratio  n.ay  be  computed  by  the  assumption  that  the  depth  of 
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soil  affected  by  the  interaction  equals  the  radius  of  the  tire. 
For  two-dimensional  conditions  the  average  strain  rate  may  be  es- 
timated as 


This  formula  was  used  to  compute  the  average  strain  rate  for 
the  test  series  on  9.00-14  tires  reported  in  Table  8 of  Ref.  1. 
For  the  tests  performed  at  the  standard  velocity  of  5 feet/second, 
the  average  strain  rate  computed  by  Eq.  (31)  varied  from  5 to 
34  percent  sec  , appreciably  higher  than  the  average  strain 
rate  in  a standard  cone  penetration  test.  For  the  tests  performed 
at  velocities  from  0.5  to  18  feet/second,  the  average  strain 
rate  computed  by  Eq.  (31)  varied  from  1 to  58  percent  sec 

Consideration  of  Strain  Rate  Effects 
in  the  Tire-Soil  Model 

In  the  test  series  on  9.00-14  tires  mentioned  previously  the 
soil  was  characteri.zed  by  its  cone  index  obtained  in  cone  penetra- 
tion tests  performed  at  the  standard  penetration  rate.  For  the 
purpose  of  model  validation,  a relationship  between  the  Coulomb 
strength  parameters  and  the  cone  index  was  established  for  the 
Buckshot  clay  in  Ref.  4,  that  was  found  to  yield  good  prediction 
for  the  tire  tests  performed  in  Buckshot  clay.  Hie  si ip -interface 
shear  parameters  nee-ded  in  the  prediction  were  established  by  an 
empirical  optimization  technique. 

In  order  to  take  strain  rate  effects  in  consideration  in  the 
model  it  would  be  desirable  to  establish  the  strength  properties 
of  the  Buckshot  clay  as  a function  of  the  strain  rate,  determine 
the  strain  rate  for  a particular  case,  and  apply  in  the  model  the 
strength  parameters  adjusted  for  the  strain  rate.  Since  the 
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strain  rate  is  not  known  a priori,  additional  iterations  would 
have  to  be  performed  in  the  model  for  its  determination.  Also, 
the  inconsistency  of  the  average  strain  rate  being  different  in  a 
standard  cone  penetration  test  from  that  in  a standard  tire  test, 
such  an  approach  would  be  unworkable.  Instead,  the  experimental 
data  reported  for  the  test  series  on  9.00-14  tires  in  Buckshot 
clay  were  analyzed  by  adjusting  the  cone  index  values  obtained  at 
the  standard  penetration  rate  for  the  variations  of  the  transla- 
tional velocity.  The  cone  penetration  resistance  was  found  to 
vary  with  the  rate  of  penetration  according  to  the  following  rela- 
tion (Ref.  30) 


Cl 

V 


0.092 


(32) 


Since  strain  rates  are  proportional  to  the  penetration  rate  in 
cone  penetration  and  to  the  translational  velocity  in  tire-soil 
interaction,  the  v/v  ratio  in  the  above  equation  may  be  substi- 
tuted  by  the  ratio  of  the  actual  translational  velocity  to  the 
standard  translational  velocity  in  tire  tests  (5  feet/second). 

An  analysis  of  tire  performance  was  made  where  the  Coulomb 
strength  parameters  were  determined  for  the  Cl  values  adjusted 
for  the  translational  velocity  by  Eq.  (32).  Results  of  this  analy- 
sis are  shown  in  Fig.  19  where  predicted  pull  coefficients  are 
shown  against  measured  ones.  Indicated  by  various  S3mbols  are  the 
ranges  of  adjustment  factors  for  the  Cl  values  used  in  the  com- 
putation. 

It  is  seen  from  Fig.  19  that  pull  coefficients  are  generally 
overpredicted  by  the  tire -soil  model  if  adjustments  in  the  strength 
properties  for  the  strain  rate  are  made  on  the  assumption  that  the 
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Fig.  19.  Comparison  of  Measured  and  Predicted  Pull  Coefficients 
for  Various  Translational  Velocities.  Predictions 
Based  on  Cone  Index  Adjusted  for  the  Ratio  of  Transla- 
tional Velocity  to  the  Standard  Velocity  (5  ft/sec)  of 
Tire  Tests 


cone  strength  varies  with  the  translational  velocities  as  indicated 
by  Eq.  (32)  . 

There  are  several  reasons  why  this  procedure  has  not  resulted 
in  a better  approximation  of  the  pull  performance.  These  can  be 
summed  up  as  follows: 

• In  the  computations  for  Fig.  19  the  interface 
friction-slip  parameters  were  assumed  on  the 
basis  of  unadjusted  Cl  values.  Th.ese  parame- 
ters may  vary  with  the  translational  velocity 
of  the  tire. 

• In  the  tire-soil  model  the  soil  properties  are 
assumed  to  be  uniform  throughout  the  affected 
soil  mass.  In  tire-soil  interaction,  strains 
and  strain  rates  vary  from  point  to  point. 

Thus,  if  variation  of  strength  properties  with 
strain  rate  is  to  be  properly  considered  in  the 
model,  then  the  model  would  have  to  be  expanded 
to  allow  for  the  variation  of  the  strength  prop- 
erties at  every  nodal  point.  This  variation 
would  have  to  be  determined  on  the  basis  of 
strain  rate  computations.  Ihe  two  methods  for 
the  computation  of  soil  inertia  forces  pre- 
sented in  Sections  IV  and  V lend  themselves 
readily  applicable  to  the  computation  of  strain 
rates.  Velocity  fields  as  well  as  particle  path 
geometries  may  be  used  for  strain  rate  computa- 
tions . 

• Hie  crude  approximation  of  the  average  strain 
rate  by  Eq.  (31)  shows  that  strain  rates  are 
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proportional  to  the  translational  velocity  for 
the  same  tire-soil  geometry  but  vary  signifi- 
cantly with  the  rut  depth  even  if  the  transla- 
tional velocity  remains  the  same.  Thus,  an 
adjustment  in  the  Cl  values  based  solely  on 
translational  velocities  cannot  reflect  the 
effect  of  strain  rate  on  the  strength  proper- 
ties in  tire-soil  interaction. 

Summary  Discussion  of  the  Effects  of  Strain  Rate 
on  Tire-Soil  Interaction  and  Conclusions 

Strain  rate  significantly  affects  the  strength  properties  of 
certain  soils  and  thereby  tire-soil  interaction.  Ihe  experimental 
information  on  the  effect  of  strain  rate  on  soil  strength  behavior 
is  limited  to  a few  t3rpes  of  soil.  The  complexity  of  the  physics 
of  strain  rate  effects  and  the  lack  of  a theory  prohibit  generaliza- 
tion of  the  results  of  experiments  for  other  types  of  soils.  Spe- 
cifically, clay  soils  saturated  to  a lesser  degree  than  in  the  ex- 
periments performed  on  Buckshot  clay  are  likely  to  exhibit  differ- 
ent behavior.  More  basic  research  is  needed  for  the  interpretation 
of  strain  rate  effects  observed  in  the  experiments. 

Theoretically,  the  most  valuable  information  on  the  effect  of 
strain  rate  on  soil  behavior  is  that  obtained  in  dynamic  triaxial 
tests  where  the  strain  rate  is  constant  and  controlled.  Indirect 
information  on  the  effects  of  strain  rate  has  been  obtained  from 
cone  penetration  and  plate  sinkage  tests  performed  at  various  rates 
of  penetration.  Interpretation  of  the  results  of  these  tests  is 
difficult  and  the  general  validity  of  correlations  with  tire  per- 
formance are  questionable  for  the  following  reasons.  First,  soil 
inertia  effects  in  these  tests  are  not  separated  from  strain  rate 
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effects.  Second,  the  strain  rate  varies  in  the  affected  soil  mass 
and,  therefore,  measured  changes  in  penetration  resistance  reflect 
the  cumulative  effect  of  strain  rate  variation  in  the  affected  soil 
volume.  Theoretical  and  experimental  research  is  needed  to  estab- 
lish strain  rate  patterns  for  these  experiments. 

Soil  strain  rates  also  vary  in  tire-soil  interaction.  To  take 
their  effect  properly  into  account  in  the  tire-soil  model  it  is 
necessary  to  determine  the  strain  rate  distribution  in  the  soil. 

Of  the  various  ways  to  accomplish  this  the  evaluation  from  parti- 
cle path  geometries  appears  to  be  the  most  practical.  More  experi- 
mental information  on  particle  path  geometries  is  needed  to  deter- 
mine the  parameters  in  the  analytical  simulation  presented  in 
Section  V. 
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VII.  USE  OF  TIRE-SOIL  MODEL  IN 
VEHICLE  RIDE  DYNAMICS  SIMULATION 


In  the  simulation  of  tire-soil  interaction  by  the  tire-soil 
model  it  is  assumed  that  tire  load  and  travel  velocity  are  constant 
and  a steady  state  exists  in  the  soil.  With  these  assumptions  the 
fundamental  relationships  governing  tire-soil  interaction  may  be 
formulated  and  compared  with  laboratory  experiments  without  undue 
interference  from  interactions  between  tire  and  vehicle. 

Another  important  area  in  mobility  research  is  the  simulation 
of  ride  d)mamics  of  off-road  vehicles . In  the  first  generation 
d5mamic  models  of  off-road  vehicles  it  is  assumed  that  the  ground 
is  rigid  and  the  motions  of  the  vehicle  are  two  dimensional.  With 
these  assumptions  it  is  possible  to  analyze  the  dynamics  of  off- 
road vehicles  under  idealized  conditions  and  validate  the  model 
for  such  conditions. 

Ideally,  in  both  vehicle  performance  and  ride  dynamics,  simu- 
lation of  the  interaction  between  running  gear  and  vehicle  should 
be  taken  into  account.  Even  though  in  the  model  expansion  reported 
in  Sections  IV  through  VI  velocity  effects  have  been  considered 
only  in  a very  preliminary  way,  the  model  may  be  used  to  demon- 
strate its  potential  use  in  vehicle  dynamic  models.  To  this  end 
sample  computations  were  performed  to  determine  axle  height-load 
relationships  for  a 9.00-14  tire  predicted  by  the  model  for  the 
conditions  shown  in  Table  5. 

Figure  20  shows  the  results  of  computations  for  Cases  1 and 
2.  The  difference  between  the  two  cases  is  in  the  value  of  the 
assumed  drawbar  pull.  In  this  case  the  100  pounds  difference 
in  drawbar  pull  does  not  appear  to  affect  the  axle  height-load 
relationship  appreciably.  However,  in  other  types  of  soils. 
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Preceding  pege  blank 


Vertical  Force  (lbs) 


TABLE  5 AXLE  HEIGHT-LOAD  RELATIONSHIPS 


Tire  Size: 
Nominal  Radius: 
Nominal  Width: 
Soil: 


9.00-14 
1.18  ft 
0.74  ft 
Buckshot  Clay 


Case  No. 

1 

2 

3 

4 

Tire  Inflation  Pressure  psi 

38 

38 

11.8 

11.8 

Limit  Pressure  psi 

29 

29 

11.6 

11.6 

Deflection  (WES  Term.) 

25% 

25% 

25% 

25% 

Deflection  Coefficient  (s) 

0.9 

0.9 

0.9 

0.9 

Translational  Velocity  ft/sec 

5 

5 

5 

10 

Drawbar  Pull  lbs 

0 

100 

0 

0 

Cone  Index  at  Standard 

37 

37 

43 

43 

Penetration  Rate 

Cone  Index  Adjusted  for 

37 

37 

43 

46 

Velocity 

especially  in  cohesionless  soils,  it  is  likely  that  the  effect  of 
the  drawbar  pull  on  the  vertical  force-axle  height  relationship 
will  be  significant.  Note  that  the  relationship  shown  in  Fig.  20 
could  be  closely  approximated  by  a linear  spring  rodel. 

Figure  21  shows  the  same  relationship  for  Cases  3 and  4 shown 
in  Table  5.  These  two  cases  differ  in  soil  strength  since  in 
Case  4 the  cone  index  was  adjusted  for  the  strain  rate  correspond- 
ing to  the  10  feet/second  translational  velocity.  The  obtained 
axle  height-load  relationships  are  almost  identical,  probably  due 
to  the  controlling  influence  of  the  limiting  pressure  that  governs 
the  normal  pressure  over  a large  portion  of  the  contact  area  in 
both  cases.  The  soil  strength  comes  into  play  only  in  a limited 
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portion  of  the  contact  area,  hence  the  almost  identical  axle 
height-load  relationships.  This  relationship  is  not  as  close  to 
a linear  one  as  the  one  shown  in  Fig.  20.  Nevertheless,  consider- 
ing all  other  idealizations  in  vehicle  d3mamic  models,  a linear 
spring  model  would  be  an  acceptable  approximation. 

These  sample  computations  show  that  the  tire-soil  model,  in 
addition  to  performance  predictions,  may  also  be  used  for  the  de- 
termination of  realistic  spring  constants  representing  both  tire 
deflection  and  ground  deformation  in  vehicle  dynamic  models.  Once 
such  constants  are  determined  for  given  tire  and  soil  parameters, 
their  inclusion  in  a dynamic  model  is  simple.  The  main  advantage 
of  this  procedure  is  that  in  the  spring  constants  determined  by 
the  tire-soil  model,  the  deformability  of  soil  is  included  allowing 
thereby  a more  realistic  simulation  of  off-road  ride  d3mamics. 
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VIII.  TOWED  PNEUMATIC  TIRE-SOIL  MODEL 

Introduction 

For  the  analysis  of  cross  country  mobility  of  all  wheel  drive 
vehicles,  determination  of  the  drawbar  pull  that  a tire  can  develop 
under  various  soil  and  terrain  conditions  is  of  prime  importance. 
Ihe  driven  pneumatic  tire-soil  model  developed  under  Contract 
DAAE07-73-C0115  serves  this  purpose.  However,  it  is  often  neces- 
sary to  evaluate  the  performance  of  off-road  vehicles  towing 
trailers  or  of  all  wheel  drive  vehicles  when  the  front  wheel  drive 
is  disengaged  or  inoperative.  For  this  purpose  it  is  necessary  to 
determine  the  towing  force  needed  to  move  a pneumatic  tire  in  vari- 
ous soils. 

For  the  determination  of  the  towing  force  a pneumatic  towed 
tire  model  was  developed  by  Janosi  in  1960  (Ref.  33)  using  soil 
parameters  from  plate  sinkage  tests.  Many  features  of  this  model 
were  incorporated  in  the  tire-soil  model  for  driven  tires.  The 
accuracy  of  performance  predictions  achieved  by  the  driven  tire 
model,  the  need  to  use  uniform  soil  parameters  in  both  driven  and 
towed  tire  models,  and  the  desire  to  incorporate  the  results  of 
more  recent  experimental  information  on  tire  geometry  and  interface 
stresses  into  the  model,  prompted  the  development  of  a towed  tire- 
soil  model  described  herein. 


Experimental  Information  on  Tov/ed 
Pneumatic  Tire  Behavior 

Numerous  towed  tire  tests  have  been  performed  by  various  re- 
searchers over  the  years.  Most  of  these  were  restricted  to  the 
measurement  of  towing  force  and  sinkage.  For  the  development  of  a 
towed  tire-soil  model  those  tire  tests  are  especially  of  interest 
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where  tire  geometry  and  interface  stresses  were  measured  (Refs.  33 
through  40) . 

The  tire  geometry  measurements  on  towed  tires  (Refs.  36  through 
40)  show  that  for  towed  tires  the  same  qualitative  relationships 
among  tire  stiffness,  soil  stiffness,  sinkage,  and  deflection  hold 
as  this  general  relationship  is  shown  schematically  in  Fig.  22  for 
driven  tires. 


Fig.  22.  Schematic  Representation  of  Tire-Soil  Behavior 
(Based  on  V/ES  Experiments) 


Experimental  information  on  stresses  measured  at  the  inter- 
face of  towed  tires  is  meager  and  restricted  to  the  measurement  of 
normal  stresses.  Since  the  salient  difference  between  driven  and 
towed  tires  is  in  the  magnitude  and  distribution  of  intarrace 
shear  stresses,  this  lack  of  experimental  information  presentv^d  a 
serious  drawback  for  the  formulation  of  interface  shear  stress 
criteria.  In  lieu  of  anything  better,  these  were  formulated  on 
the  basis  of  interface  shear  stress  measurements  on  towed  rigid 
wheels.  Figure  23  shows  a typical  result  of  such  measurements  by 
Krick  (Ref.  40) . 


Pig.  23.  Distribution  of  Normal  and  Shear  Stresses  Beneath 
a 0.88  in.  Diameter  Rigid  V/heel  (Averaged  Over  the 
V/idth  for  Two  Tests)  Wheel  Load:  6i^0  kg  (Prom  Ref. 


Concept  c2  Towed  Tire-Soil  Interaction 


In  the  conceptual  development  of  the  towed  tire-soil  model 
the  inferences  drawn  from  the  study  of  experimental  information 
were  incorporated  in  the  following  way. 

• Tire  deflection  was  related  to  the  arc 
length  of  the  front  slip  line  field  in 
the  same  v7ay  as  in  the  driven  tire-soil 
model. 

• The  interface  friction  coefficient,  6,  that 

governs  the  shear  stress  distribution,  was 
assumed  to  vary  linearly  from  the  entry  and 
exit  points  to  a central  angle,  where 

its  sign  changes.  This  assumption  results  in 
a shear  stress  distribution  similar  to  that 
obtained  for  rigid  wheels  and  shown  in  Fig.  23. 

• For  the  shear  stresses  the  criterion  was 
adopted  that  their  moment  about  the  tire  axis  be 
zero.  While  at  first  glance  this  criterion  ap- 
pears to  be  self-evident  for  a.  towed  tire  where 
the  applied  torque  is  zero,  it  should  be  brought 
to  mind  that  in  the  case  of  a deflecting  tire 
normal  stresses  may  also  produce  moments  about 
the  axis.  In  such  a case  the  moment  produced 
by  the  normal  stresses  would  have  to  be  balanced 
by  the  shear  stresses  resulting  in  a different 
shear  stress  distribution.  It  is  in  this  respect 
that  the  lack  of  experimental  information  on  in- 
terface shear  stresses  beneath  towed  tires  hin- 
dered the  development  of  a towed  tire  model  that 
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would  be  responsive  to  this  contingency  in 
shear  stress  distribution. 

• For  the  driven  tire-soil  model,  relationships 
were  developed  between  slip  and  the  interface 
friction  angle.  For  towed  tires,  skid  or 
negative  slip  is  of  less  importance  since  it 
does  not  affect  travel  efficiency  or  power 
consumption.  Therefore,  skid-interface  fric- 
tion relationships  were  not  analyzed  and 
included  in  the  model. 

Development  of  Towed  Tire-Soil  Model 

The  concepts  of  towed  tire-soil  interaction  described  previ- 
ously were  incorporated  in  the  towed  tire -soil  model  shown  in 
Fig.  24.  The  tire  centerline  geometry  consists  of  logarithmic 
spirals  and  straight  lines,  as  in  the  driven  tire  model.  Tire 
deflection  is  assumed  to  start  ahead  of  the  entry  point  and 
ceases  after  the  exit  point.  This  assumption  is  incorporated  in 
the  model  by  the  angle  a'  (see  Fig.  24)  that  defines  the  end 
points  of  the  logarithmic  spirals  in  the  front  and  rear.  The 
normal  stresses  are  controlled  by  plastic  equilibrium  c'jnditions 
in  the  soil  whenever  they  are  less  than  the  limit  pressure  p^ 
that  depends  on  the  tire  inflation  pressure.  The  shear  stresses 
and  the  interface  friction  angle  6 are  assumed  to  vanish  at  the 
central  angle  a^.  A linear  variation  of  6 from  the  entry  and 
exit  points  results  Ln  a shear  stress  distribution  as  shown  in 
Fig.  24. 

For  towed  tire-soil  interaction,  the  soil  model  is  the  same 
as  for  driven  tire-soil  interaction.  Soil  behavior  is  represented 
by  its  Coulomb  strength  parameters  v'lnd  its  unit  weight.  The  gen- 
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Pig.  2^.  Towed  Tire-Soil  Model 

eral  concept  of  towed  tire-soil  interaction  as  incorporated  in  the 
towed  tire-soil  model  shown  in  Fig.  24  does  not,  in  itself,  uniquely 
define  the  geometry  of  the  tire-soil  interface.  The  condition  that 
the  normal  stress,  q at  the  assumed  angle,  o.  , be  equal  to 

the  limit  pressure  p^  defines  the  entry  angle,  since  it  is  by 
the  variation  of  the  entry  angle  that  this  condition  is  met.  How- 
ever, the  rear  angle  is  undefined.  Experimental  information 

on  u^e  geometry  of  towed  tiic«=  shows  that  the  rear  angle  varies 
with  the  load  and  the  entry  angle.  To  define  a general  relationship 


for  a wide  variety  of  soil  and  loading  conditions  the  sasne  approach 
was  followed  as  for  the  development  of  the  driven  tire-soil  model. 

A computer  program  was  written  for  the  towed  tire-soil  model  and 
towing  force  computations  were  made  for  a variety  of  assumptions 
regarding  the  rear  angle.  The  towed  force  computations  were  per- 
formed for  tire  sizes  and  soil  conditions  for  which  experimental 
data  were  reported  in  Ref.  1.  These  towed  tire  tests  were  selected 
for  validation  of  the  model  because  soil  strength  parameters  for 
the  two  t3T3es  of  soils  used  in  the  tests  were  already  established 
in  Report  No.  11900  (LL147)  (Ref.  4) . 

Since  the  cost  of  the  large  number  of  computer  runs  required 
to  perform  the  contemplated  systematic  analysis  of  the  effect  of 
the  various  assumptions  on  the  towed  force  is  high,  the  computer 
program  was  adapted  to  the  General  Data  Corporation  Nova  Minicom- 
puter owned  by  the  Research  Department.  Running  time  of  the  com- 
puter program  on  the  minicomputer  is  about  4-8  times  that  on  large 
computers  but  the  running  cost  is  nominal.  Arrangements  were  made 
for  input  and  output  data  readiug  from,  and  recording  on,  data 
files  so  that  results  for  a large  number  of  selected  cases  could 
be  obtained  overnight. 

The  computer  analyses  of  a large  number  of  cases  served  the 
dual  purpose  of  comparing  results  with  experimental  data  and  check- 
ing on  tolerance  limits  and  convergency  criteria  of  various  itera- 
tion schemes  adopted  in  the  computer  program.  In  connection  with 
the  latter,  some  new  problems  arose  that  are  discussed  below. 

Problems  ISncountered 

In  the  development  of  the  computer  program  for  towed  tire -soil 
interaction  the  experiences  gained  with  the  driven  tire-soil  model 
were  heavily  relied  upon.  One  of  the  tasks  that  had  to  be  solved 
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in  both  model?  is  that  of  finding  the  size  of  the  slip  line  field 
that  matches  a prescribed  normal  stress  at  the  interface  at  its 
end  point.  All  iteration  schemes  developed  for  the  driven  tire- 
soil  model  were  based  on  the  assumption  that  the  normal  stresses 
computed  for  any  slip  line  field  increase  from  the  edge  point 
adjoining  the  free  surface  toward  the  inner  points.  While  this  is 
generally  true  for  a constant  6 angle  of  interface  friction 
assumed  for  driven  tires,  it  is  not  necessarily  so  when  the  angle 
decreases  from  the  edge  toward  the  inner  point.  It  was  found  that 
in  certain  cases  in  cohesive  soils  with  a low  friction  angle  the 
normal  stress  is  unchanged  or  even  decreases  from  the  edge  point 
for  a certain  rate  of  decrease  in  the  6 angle.  Measurements  of 
normal  stresses  at  the  interface  of  a rigid  wheel  towed  in  clay 
reported  by  Uffelmann  (Ref.  41)  and  shown  in  Fig.  25  indicate  that 
in  highly  cohesive  soils  this  is,  indeed,  possible. 


Pig.  25.  Distribution  of  Normal  Stresses  Beneath  a 5^-ih.  Diameter 
Rigid  V/heel  Towed  in  Heavy  Clay.  Load:  6720  lbs;  Slip; 

= 27^,  (From  Ref.  4l) 
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This  feature  of  the  normal  stress  distribution  in  cohesive 
soils  made  it  necessary  to  make  provisions  in  the  program  that 
account  for  this  situation.  The  main  problem  that  had  to  be  re- 
solved was  the  determination  of  the  size  of  the  front  slip  line 
field  for  nearly  constant  normal  stresses.  In  principle,  this 
problem  was  resolved  by  assuming  that  the  size  of  the  front  field 
for  entry  angles  between  30  and  60  degrees  is  determined  by 
the  condition  that  the  vertical  components  of  all  stresses  be 
equal  to  the  load.  The  rear  angle  is  assumed  as  10  degrees  for 
the  variations  of  the  entry  angle  from  30  to  60  degrees.  For 
entry  angles  less  than  30  degrees  the  rear  angle  is  reduced  by 
half  of  the  difference  between  the  entry  angle  and  30  degrees.  A 
computational  scheme  corresponding  to  these  assumptions  was  incor- 
porated in  the  computer  program.  Results  of  the  towing  force  com- 
putations were  compared  with  experimental  data  and  were  found  to 
yield  acceptable  approximations  within  the  limitations  of  a two 
dimensional  model. 

Results  of  the  Analyses  of  Experimental  Data 

The  systematic  analyses  of  the  experimental  data  were  directed 
primarily  toward  determination  of  the  best  assumption  as  to  the 
relative  positions  of  the  front  and  rear  slip  line  field  controlled 
by  the  angle  of  separation  and  rear  angle.  The  effects  of  other 
assumptions,  such  as  the  relationship  between  the  limiting  and  in- 
flation pressures,  the  value  of  the  interface  friction  angle  at  the 
entry  and  exit  point,  and  the  location  of  zero  shear  stress  were 
also  investigated.  In  the  course  of  these  systematic  analyses  per- 
formed on  the  minicomputer  over  3000  cases  were  investigated,  all 
representing  different  combinations  of  the  input  data  and  assump- 
tions. The  cost  of  this  systematic  investigation,  if  carried  out 
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on  large  computers  at  commercial  rates,  would  have  amounted  to 
more  than  half  of  the  contract  cost.  Even  though  the  number  of 
cases  investigated  appears  to  be  large,  only  a limited  number  of 
all  possible  combinations  of  input  variables  and  assumptions  are 
covered.  Nevertheless,  it  was  possible  to  establish  approximate 
relationships  that  define  the  towed  tire-soil  model  completely  and 
result  in  acceptable  prediction  accuracy.  Although  further  sys- 
tematic analyses  could  have  resulted  in  refinements  in  the  model, 
the  time  spent  would  have  been  prohibitive.  Also,  the  systematic 
analyses  showed  that  prediction  accuracy  could  not  be  improved  sig- 
nificantly unless  three  dimensional  conditions  were  considered  in 
the  model. 

The  following  assumptions  were  incorporated  in  the  final  com- 
puter program  on  the  basis  of  these  systematic  analyses: 

Initial  value  of  rear  angle  ” 10® 

Initial  value  of  separation  angle  = 20® 

Relationship  between  rear  and  separation  angle  = 

a = 1.5  X a +5°  (33) 

In  IT 

The  effects  of  limiting  pressure  on  the  prediction  of  towing 
forces  were  also  investigated.  In  the  case  of  towed  tires  the 
following  relationship  was  found  to  yield  good  simulations 


p(psi)  = Q^82  X p,  + 6 


The  interface  friction  angle  6 is  assumed  to  vary  from  a 
maximum  value  at  the  entry  and  exit  points  to  0 at  some  inner 
point.  The  variation  of  6 along  the  interface  is  thus  defined 
by  its  maximum  value  and  the  angle  at  which  its  value  becomes 

zero.  This  angle  was  defined  in  the  program  as 
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a = F a 
o T m 


(35) 


where  the  value  of  the  factor  was  to  be  determined  by  the  con- 

dition that  the  moment  of  the  shear  stresses  about  the  tire  axis 
be  zero.  The  systematic  analyses  were  run  with  various  values  of 
the  factor.  It  was  found  that  a value  of  F^  = 0.5  generally 
yielded  a negligible  moment  of  the  shear  stresses.  Although  it 
would  have  been  possible  to  set  up  an  iteration  scheme  whereby  the 
moment  of  the  shear  stresses  was  made  less  than  a specified  tol- 
erance by  an  appropriate  variation  of  the  F^  value,  it  was  found 
that  the  prediction  accuracy  would  not  have  been  improved  in  a sig- 
nificant degree  to  warrant  the  additional  computer  time  required 
for  these  iterations.  If,  however,  for  any  reason,  the  moment  of 
shear  stresses  appears  to  be  unacceptably  high,  an  appropriate 
change  in  the  input  F^  value  would  bring  it  to  an  acceptable 
level. 

The  maximum  value  of  the  interface  friction  angle  was  assumed 
as 


6 

max 


(36) 


In  the  systematic  analysis  the  effect  of  the  factor  Fg  on  the 
towed  force  coefficient  was  investigated  by  assuming  various  values 
for  Fg.  It  was  found  that  for  good  simulation  Fg  may  vary  be- 
tween 0.25  and  0.50.  These  values  are  in  the  range  of  Fg 
values  evaluated  from  experimental  data  on  rigid  wheels.  Shear 
stress  measurements  on  towed  pneumatic  tires  in  various  soils  would 
be  needed  to  firm  up  the  value  of  the  Fg  factor.  A detailed  flow 
chart  of  the  computer  program  incorporating  these  features  is  given 
in  Appendix  C. 
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A comparison  of  measured  and  predicted  towed  force  coeffi- 
cients is  given  in  Tables  6 through  10  for  various  tire  sizes, 
loading,  and  soil  conditions.  The  experimental  data  were  reported 
in  Ref.  1 where  more  detailed  information  on  the  performance  of 
experiments  may  be  found. 

As  shown  in  the  tabulations,  in  addition  to  the  towed  force 
coefficient  actually  measured  in  the  experiment,  the  one  computed 
from  the  best  fitting  curve  for  all  experiments  performed  in  the 
same  soil,  as  reported  in  Ref.  1,  is  also  given.  Since  all  ex- 
perimental measurements  necessarily  contain  some  experimental  error, 
differences  in  predicted  and  measured  towed  coefficients  may  arise 
from  errors  or  inaccuracies  in  either  value.  In  the  course  of  the 
systematic  analyses  of  the  experimental  data  it  became  evident,  in 
some  instances,  that  efforts  to  simulate  individual  experiments 
may  be  antiproductive  if  that  experiment  contains  a significant 
error.  A comparison  of  measured  towed  coefficients  with  that  ob- 
tained from  a best  fitting  curve  is  an  indication  of  whether  that 
may  be  the  case. 

The  predicted  towed  force  coefficients  were  obtained  with 
F =0.5  and  F,.  = 0.25  and  deflection  coefficients  e estab- 

T 6 

lished  for  driven  tires  in  Ref.  4 and  shown  in  Table  11. 


Figures  26  and  27  show  typical  centerline  geometries  and  in- 
terface stress  distributions  predicted  by  the  model  for  the  follow- 
ing conditions: 


Tire  Size; 

Tire  radius: 

Tire  Width: 

Inflation  Pressure: 
Deflection  Coefficient 
Soil: 

Cone  Index: 

Cone  Index  Gradient: 


e: 


Fig.  26 

9.00-14 
1.18  ft 
0.73  ft 
12.0  psi 
0.905 
Yuma  Sand 

3.5  pci 


Fig.  27 

4.00-20 
1.17  ft 
0.37  ft 
40.0  psi 
0.945 

Buckshot  Clay 
18  psi 
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TABLE  6 COMPARISON  OF  PREDICTED  AND  >EASL’RED  TOWED  FORCE  COEFFICIENTS 

TIRE  SIZE:  6,00-16 

SOIL:  YUMA  SAND 


Case 

Code 

No. 

WES 

Idenciflcaclon 

No. 

Cone  Index  (psl) 
(In  clay)  or 
Cone  Index 
Gradient  (pci) 
(In  sand) 

Load 

lbs 

Inflation 

Pressure 

psl 

WES 

Deflection 

6/h 

Towed 

Force  Coefficient 

Measured 

Averaged 
From  Best 
Fitting  Curve 

Predicted 

I 

164 

802A 

5.8 

213 

8.5 

0. 

15 

0.155 

0.114 

0.088 

2 

164 

805A 

U.l 

215 

8.5 

0.070 

0.058 

0.081 

3 

164 

809A 

15.6 

222 

8.5 

0.059 

0.048 

0.075 

4 

164 

808A 

12.4 

293 

11.4 

0.065 

0.073 

0.085 

5 

164 

807A 

10.4 

458 

17.2 

0.131 

0.139 

0.121 

6 

165 

35A 

4.0 

650 

29.0 

0.449 

2.3 

0.203 

7 

164 

816A 

13.5 

240 

4.5 

0. 

25 

0.042 

0.0395 

0.078 

8 

165 

37A 

15.0 

223 

4.5 

0.063 

0.035 

0.083 

9 

164 

818A 

15.6 

455 

10.3 

0.040 

0.057 

0.092 

mm 

165 

33A 

2.3 

429 

10.3 

0.424 

0.604 

0.238 

11 

164 

812A 

13.8 

865 

21.0 

0.088 

0.116 

0.127 

12 

164 

817A 

9.5 

863 

21.0 

0.200 

0.179 

0.155 

13 

164 

803A 

5.8 

225 

2.5 

0 

35 

>>.116 

0.055 

0.111 

14 

164 

813A 

16.4 

239 

2.5 

0.046 

0.078 

15 

164 

8UA 

17,3 

446 

7,0 

0.018 

0.092 

16 

165 

34A 

3.5 

674 

10.3 

0,365 

0.217 

17 

164 

Alii 

15.0 

870 

13.0 

0.055 

0.125 

SOIL: 

BUCKSHOT  CUY 

32 1C 

20 

238 

8.5 

0. 

15 

0.059 

0.087 

iSI 

323C 

63 

237 

8.5 

0,013 

H.S.F. 

wm 

329C 

40 

238 

8.5 

0.021 

H.S.F. 

167 

32  5C 

20 

430 

17.2 

0.144 

0.142 

0.124 

168 

351C 

55 

461 

17.2 

0,033 

0.042 

H.S.F. 

169 

327C 

20 

631 

29.0 

0,254 

0.367 

0.157 

170 

339C 

40 

688 

29.0 

0.112 

0.093 

0.086 

mm 

344C 

37 

877 

38.0 

J 

0.155 

0.170 

0.125 

172 

359C 

54 

875 

38.0 

f 

0.114 

0.089 

H.S.F. 

173 

322C 

20 

243 

4.5 

25 

0.037 

0.047 

H.S.F. 

174 

324C 

61 

234 

4.5 

0.004 

0.097 

H.S.F. 

175 

330C 

42 

237 

4.5 

0.021 

0.051 

H.S.F. 

176 

326C 

19  ' 

44>i 

10.3 

0.119 

0,107 

0.135 

177 

357C 

51 

460 

10.3 

0.022 

0,036 

H.S.F. 

178 

354C 

50 

662 

15.3 

0.050 

0.051 

H.S.F, 

179 

338C 

22 

688 

15.3 

0.222 

0.179 

0.165 

180 

362C 

50 

716 

16.2 

0.052 

0.056 

H.S.F. 

181 

337C 

21 

848 

21.0 

0.292 

0.355 

0.183 

182 

345C 

37 

893 

21.0 

0.127 

0.111 

H.S.F. 

183 

360C 

52 

889 

21.0 

0.066 

0,069 

H.S.F. 

184 

33 1C 

20 

240 

2.5 

Bl 

35 

0.0::8 

0.040 

H.S.F. 

185 

342C 

37 

226 

2.5 

0.027 

0.024 

H.S.F. 

186 

355C 

53 

223 

2.5 

0.045 

0.019 

H.S.F. 

187 

332C 

20 

457 

7.0 

0.074 

0.081 

0.092 

Big 

343C 

37 

448 

7.0 

0.033 

0.040 

H.S.F. 

ICTl 

356C 

53 

447 

7.0 

0.034 

0.0295 

H.S.F. 

190 

34 1C 

36 

672 

■SSh 

0.042 

0.063 

H.S.F. 

191 

338C 

50 

673 

BSbH 

0.022 

0.044 

H.S.F. 

197 

3 oC 

39 

873 

0.071 

0.079 

H.S.F. 

193 

361C 

52 

87) 

0.046 

0.056 

H.S.F. 

■ — - 

— 1 

1 

TAfo'^  7 C»r.v'ia.SON  OF  PREDICTED  ASt)  MEASliRKD  TOWD  FORCE  COEFFICIENTS 

TIKE  SIZE;  9.00-14 


son.;  YUMA  SAND 


WES 

Idonclflcaclon 

No. 


164  778A 
164  779A 

164  780A 
164  786A 
164  777A 
164  782A 
164  783A 
164  785A 

164  781A 
164  784A 


Cone  Indi>x  (psl) 
(In  clay)  or 
Cone  Index 
Gradient  (pci) 
(In  sand) 


Inn.iclon 

Pressure 

psl 


WES 

Deflection 

6/h 


Towed  Force  Coefficient 


■Measured 


Averaged 
From  Best 
Fitting  Curve 


165  28A 

165  3 A 

165  24A 

1-65-64 
1-65-65 
1-65-66 
1-65-67 
1-65-68 
1-65-69 
1-65-70' 
1-65-71 
1-65-72 
1-65-74 

165  9A 
165  llA 


TABLE  8 COMPARISON  OF  PREDICTED  AND  HEASlfRED  TOWED  FORCE  COEFFICIENTS 

TIRE  SIZE:  4.00-7 
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SOIL:  YIU\  SASD 


Case 

Code 

No. 

WES 

Idenclflcatlon 

No. 

Cone  Indo.x  (psl) 
(In  clay)  or 
Cone  Index 
Gradient  (pci) 
(In  sand) 

Load 

lbs 

Int latlon 
Pressure 
psl 

WES 

Deflection 

'/h 

Tin  . .! 

•*.»ivo  Co«*t  i li’ 

io.it 

Heasurcc 

Averae.ed 
From  Best 
Fitting  Curvt 

Predicted 

44 

164 

198A 

17.3 

85 

16.2 

0 

.15 

0.259 

0.085 

0.068 

45 

164 

824A 

15.8 

106 

16.2 

0.113 

0.117 

0.082 

46 

164 

82  5 A 

9.8 

123 

33.2 

0.171 

0.258 

0.250 

47 

164 

800A 

13.0 

210 

33.2 

0.305 

0.377 

0.175 

48 

164 

827A 

17.6 

121 

6.2 

0 

.25 

0.025 

0.073 

0.107 

49 

164 

828A 

19.6 

122 

6.2 

0.033 

0.067 

0.104 

50 

164 

831A 

25.1 

185 

17.0 

0.130 

0.077 

0.098 

51 

164 

822A 

13.0 

216 

17.0 

0.139 

0.188 

0.161 

52 

164 

829A 

19.6 

234 

17.0 

0.107 

0.126 

0.131 

53 

164 

826A 

12.1 

348 

26.0 

0.210 

0.434 

0.229 

54 

164 

833a 

18.7 

109 

2.5 

0 

.35 

0.119 

0.049 

0.100 

55 

164 

834a 

18.7 

152 

5.5 

0.086 

0.064 

0.109 

56 

164 

lA 

22.8 

145 

5.5 

0.103 

0.052 

0.100 

57 

164 

830A 

19.6 

224 

10. 1 

0.076 

0.085 

0.132 

58 

164 

832a 

19.9 

455 

12.7 

r 

0.125 

0.174 

0.157 

SOIL:  BUCKSHOT  CLAY 

21 

363C 

56 

113 

16.2 

0 

.15 

0.027 

0.034 

H.S.F. 

22 

36  8C 

45 

117 

16.2 

0.117 

0.040 

H.S.F. 

23 

371C 

26 

97 

16.2 

0.093 

0.062 

0.079 

24 

365C 

46 

226 

33.2 

0.137 

0.087 

0.077 

25 

383C 

66 

228 

33.2 

0.083 

0.056 

H.S.F. 

26 

373C 

41 

335 

51.8 

0.188 

0.212 

0.125 

27 

380C 

66 

340 

51.8 

0.109 

0.092 

H.S. 

28 

375C 

42 

439 

63.5 

0.237 

0.427 

0.143 

29 

385C 

66 

451 

63.5 

0.146 

0.141 

H.S.F. 

30 

364C 

46 

117 

6.2 

0.25 

0.026 

0.033 

H.S.F. 

31 

372C 

26 

103 

6.2 

0.068 

0.050 

H.S.F. 

32 

36  7C 

62 

217 

17.0 

0.134 

0.044 

H.S.F. 

33 

370C 

44 

229 

17.0 

0.074 

0.068 

H.S.F. 

34 

374C 

42 

330 

26.0 

0.148 

0.122 

H.S.F. 

35 

384  C 

66 

340 

26.0 

0.097 

0.067 

H.S.F. 

36 

376c 

37 

436 

35.0 

0.287 

0.292 

0.166 

386c 

65 

446 

35.0 

0.108 

0.097 

H.S.F. 

38 

378C 

34 

221 

10.4 

0 

.35 

0.104 

0.073 

H.S.F. 

39 

387C 

66 

222 

10.4 

0.072 

0.037 

H.S.F. 

40 

379C 

36 

322 

17.0 

0.203 

0.115 

H.S.F. 

41 

388C 

66 

336 

17.0 

0.042 

0.055 

H.S.F. 

42 

377C 

38 

446 

21.9 

0.197 

0.192 

H.S.F. 

43 

389C 

68 

448 

21.9 

0.083 

0.074 

H.S.F. 
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TABLE  9 COMPARISON  OF  PREDICTED  AND  MEASURED  TOWED  FORCE  COEFFICIENTS 

TIRE  SIZE:  4.00-20 


SOIL:  YUMA  SAND 


Case 

Code 

No. 

WES 

Identification 

No. 

Cone  Index  (psl) 
(in  clay)  or 
Cone  Index 
Gradient  (pci) 
(in  sand) 

1 

Inflation 

Pressure 

psi 

WES 

Deflection 

6/h 

Towed 

Force  Coefficient 

Measured 

Averaged 
From  Best 
Fitting  Curve 

Predicted 

59 

164 

791A 

7.5 

218 

24.7 

0.15 

0.211 

■8H 

60 

164 

793A 

16.4 

221 

24.7 

0.118 

61 

164 

738a 

12.1 

426 

48.2 

0.246 

0.119 

62 

164 

794  A 

12.7 

440 

48.2 

0.227 

0.109 

63 

164 

795A 

16.1 

446 

48.2 

0.195 

0.181 

0.093 

64 

165 

14A 

27.7 

22? 

11.4 

0 

.25 

0.079 

0.038 

0.076 

65 

165 

15A 

14.4 

238 

11.4 

0.046 

0.065 

0.091 

66 

165 

19A 

16.1 

338 

18.2 

0.062 

0.079 

0.100 

67 

165 

16a 

13.0 

450 

24.7 

i 

0.156 

0.131 

0.123 

68 

165 

21A 

25.6 

233 

6.7 

0 

.35 

0.064 

0.034 

0.084 

69 

165 

22A 

25.6 

350 

11.0 

0.029 

0.058 

0.0*11 

70 

16  5 

20a 

16.4 

442 

15.0 

1 

0.048 

0.074 

O.lll 

SOIL 

BUCKSHOT 

CLAY 

105 

401C 

48 

305 

82.0 

0 

tm 

0.056 

0.065 

H.S.F. 

106 

402C 

22 

307 

82.0 

0 

.08 

0.156 

0.251 

0.123 

107 

26  9C 

20 

204 

24.7 

0 

.15 

0.059 

0.086 

0.099 

108 

27SC 

48 

221 

24.7 

0.023 

0.036 

H.S.F. 

109 

279c 

32 

228 

24.7 

0.048 

0.055 

0.067 

110 

271C 

18 

388 

40.0 

0.255 

0.427 

0.177 

111 

273C 

46 

456 

48.2 

0.070 

0.081 

0.085 

112 

287C 

32 

448 

48.2 

0.127 

0.139 

0.125 

113 

286c 

48 

629 

61.0 

0.097 

0.122 

0.099 

114 

270C 

22 

204 

11.4 

0.25 

0.044 

0.057  ■ 

0.099 

115 

276c 

52 

228 

11.4 

0.028 

H.S.F. 

116 

280C 

33 

225 

11.4 

0.042 

H.S.F. 

117 

272c 

19 

377 

18.2 

KSm 

0.171 

0.162 

118 

274C 

48 

450 

24.7 

0.056 

H.S.F. 

119 

288C 

33 

447 

24.7 

0.103 

0.114 

120 

281C 

47 

645 

37.5 

0.095 

121 

291C 

40 

664 

37.5 

0.130 

122 

293C 

54 

637 

37.5 

0.055 

0.073 

H.S.F. 

123 

277C 

50 

211 

6.7 

■ 

.35 

0.03E 

0.025 

H.S.F. 

124 

283C 

29 

236 

6.7 

0.021 

0.042 

H.S.F. 

125 

289C 

19 

236 

6.7 

0.021 

0.067 

0.110 

126 

278c 

50 

433 

15.0 

0.035 

0.044 

H.S.F. 

127 

284C 

30 

454 

15.0 

0.077 

0.125 

128 

295C 

19 

421 

15.0 

0.140 

0.181 

129 

282C 

45 

651 

23.0 

0.060 

0.078 

H.S.F. 

130 

292C 

35 

664 

23.0 

0.092 

0.115 

0.136 

131 

294C 

52 

628 

23.0 

0.062 

H.S.F. 

132 

285C 

52 

696 

23.0 

0.052 

0.070 

H.S.F. 

133 

400C 

48 

660 

16.1 

0.45 

0.055 

0.063 

H.S.F. 

134 

403C 

22 

635 

16.1 

- 

0.228 

0.210 

0.183 
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Fig.  26.  Interface  Normal  and  Shear  Stress  in  Sand  Predicted 
by  the  Towed  Tire-Soil  Model 

Evaluation  of  the  Towed  Tire-Soil  ^fodel 

In  Tables  6 through  10  towed  force  coefficients  predicted  by 
the  model  are  compared  with  experimental  data.  The  soil  strength 
parameters  used  in  these  predictions  are  those  established  for 
Yuma  StiUd  and  Buckshot  clay  in  Ref.  4 as  functions  of  cone  penetra- 
tion values.  Considering  that  the  soil  strength  parameters  deter- 
mined on  this  basis  are  approximate  ones  the  prediction  accuracy 
achieved  by  the  towed  tire-soil  model  is  generally  very  good. 

A review  of  towed  force  coefficient  predictions  was  made  with 
the  objective  of  determining  causes  of  differences  between  predic- 
tions and  experiments.  As  a result  of  this  review  the  following 
conclusions  were  reached. 


Pig,  27.  Interface  Normal  and  Shear  Stresses  in  Clay  Predicted 
by  the  Towed  Tire-Soil  Model 

• Underpredictions  of  the  towed  force  coefficient 
occurred  when  the  measured  coefficient  was  very 
high  (higher  than  about  0.3).  The  major  cause 
of  this  underprediction  was  that  three  dimen- 
sional effects  were  disregarded  in  the  two  dimen- 
sional model.  Soil  failure  in  the  third  (lateral) 
dimension  would  limit  normal  stresses  in  the 
transverse  direction  beneath  the  tire  and  result 
in  a reduction  of  normal  stresses  toward  the  side 
edges  of  the  tire,  especially  in  the  front  field. 
Because  of  this  reduction  of  the  normal  stresses 
the  contact  area  would  have  to  increase  to  bal- 
ance the  load  resulting  in  higher  entry  angle, 
sinkage  and  towed  force  coefficient. 
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• Bow  waves,  unaccounted  for  in  the  model,  may 
also  cause  an  increase  in  the  towed  force 
coefficient. 

• In  cohesive  soils,  the  towed  force  coeffi- 
cient is  very  sensitive  to  small  changes  in 
soil  strength.  This  may  be  best  illustrated 
by  the  summary  plot  of  the  WES  tire  tests 
(Fig.  28)  where  test  results  were  plotted 
against  the  dimensionless  number 


N 


c 


Cl  • 


b » d /6\2 
W Vh/ 


1 + b/2d  • 


When  N is  less  than  about  three,  the  cowed 
c 

force  coefficient  curve  becomes  very  steep 
indicating  that  very  small  changes  in  the 
cone  index  result  in  large  changes  in  the 
towed  force  coefficient.  Under  these  condi- 
tions it  would  be  unrealistic  to  expect  that 
the  towed  force  coefficient  could  be  deter- 
mined by  any  method  very  accurately. 

• The  towed  tire-soil  model  is  based  on  the 
concept  that  soil  failure  conditions  govern 
the  towed  tire-soil  interaction.  In  strong 
cohesive  soils,  soil  failure  does  not  occur 
under  the  tire  load  and  the  concept  of  tire- 
soil  interaction  does  not  apply.  In  the 
computer  program  this  situation  is  recognized 
by  an  approximate  calculation  that  computes 
the  failure  stresses  for  an  infinitesimally 
small  slip  line  field  and  compares  these 
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Towed  Force  Coefficient 


Fig,  28.  Variation  of  Tovied  Force  Coefficients  with  in 
Fat  Clay  (From  Ref.  1) 


with  the  limit  stress.  If  these  computed 
failure  stresses  are  larger  than  the  limit 
stress  then  no  failure  occurs  in  the  soil, 
and  therefore,  soil  failure  is  not  the  con- 
trolling factor  in  tire-soil  interaction. 
The  computer  message  for  this  condition  is 
"use  hard  surface  formula"  indicated  by 
"H.S.F."  in  the  tabulation. 


In  the  towed  tire  tests  on  Buckshot  clay, 
hard  surface  condition  is  frequently  indicated 
in  the  predicted  towed  force  coefficient  column. 
An  examination  of  these  cases  showed  that  in 
the  majority  of  the  cases  the  measured  towed 
force  coefficient  was  in  the  range  of  0.02-0.05. 
Since  the  rolling  resistance  of  tires  on  con- 
crete or  asphaltic  surface  is  in  the  lower  limit 
of  this  range,  the  hard  surface  designation 
appears  to  be  appropriate.  The  measured  sinkage 
in  these  cases  is  very  small  indicating  that  the 
soil  actually  behaved  as  a stiff  rather  than  a 
plastic  material. 

Tire-soil  interaction  in  these  cases  is 
goverried  by  the  deformation  rather  than  by  the 
strength  properties  of  soil.  For  the  simula- 
tion of  tire-soil  interaction  under  such  condi- 
tions, the  soil  model  would  have  to  be  modified 
accordingly.  Since  this  situation  is  not  criti- 
cal in  off-road  mobility  the  development  of  such 
a model  was  not  a tempted. 


IX.  CONCLUSIONS  AND  RECOMMENDATIONS 


Methods  have  been  developed  to  estimate  soil  inertia  forces 
generated  by  a tire  traveling  at  various  velocities  and  to  take 
these  inertia  fo-  :es  into  account  in  the  solutions  of  the  differen- 
tial equations  of  plasticity  for  soils.  These  methods  have  been 
incorporated  in  the  computer  program  simulating  tire-soli  interac- 
tion. From  the  analysis  of  the  effect  of  soil  inertia  forces  on 
tire-soil  interaction  the  following  conclusions  were  drawn. 

As  long  as  the  inertia  forces  generated  in  the  soil  by  the 
traveling  tire  result  in  single-valued  solutions  for  the  front  and 
rear  slip  line  field,  the  effect  of  soil  inertia  forces  on  tire- 
soil  interaction  is  minor.  At  a relatively  low  velocity  (in  the 
range  of  10  feet /second)  solution  of  the  differential  equations 
of  plasticity  becomes  multivalued  indicating  that  soil  behavior 
must  be  different  from  that  described  by  plasticity  theory.  The 
effect  of  such  stress  states  on  soil  behavior  is  unexplored  and 
cannot  be  analyzed  by  present  state  of  the  art  methods. 

It  is  recommended  that  both  theoretical  and  experimental  re- 
search be  started  to  gain  i.nsight  in  soil  behavior  under  such  con- 
ditions and  a theory  be  developed  that  could  be  applied  to  tire- 
soil  interaction  at  high  speeds. 

In  tire-soil  interaction  the  rate  of  loading  of  the  soil  and 
the  rate  of  strain  is  directly  affected  by  the  velocity  of  travel. 
The  strength  properties  of  certain  soils  are  appreciably  affected 
by  changes  in  the  rate  of  strain,  while  other  soils  are  relatively 
insensitive  to  strain  rate  variations.  The  increase  of  soil 
strength  with  strain  rate  is  probably  the  major  contributor  to  the 
improvement  of  tire  performance  with  speed. 
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strain  rates  vary  in  the  affected  soil  mass  in  both  tire- 
soil  interaction  and  in  the  various  field  tests  used  for  soil 
property  determination.  In  order  to  approach  the  problem  of 
strain  rate  effects  in  a logical  way,  it  is  necessary  to  develop 
methods  for  the  determination  of  strain  rates  in  tire-soil  inter- 
action as  well  as  in  the  various  field  tests. 

Strain  rates  are  directly  related  to  particle  path  geometry. 

It  is  recommended  that  methods  be  developed  for  the  determination 
of  strain  rates  from  particle  paths  and  experimental  determination 
of  particle  path  geometries  be  extended  to  tire-soil  interaction, 
cone  penetration,  and  plate  sinkage  tests  under  a variety  of  soil 
conditions.  More  experimental  information  on  particle  path  geome- 
tries would  also  be  useful  to  improve  the  computation  of  soil 
inertia  forces  in  tire-soil  interaction. 

An  important  new  concept,  that  of  the  contact  slip,  was  de- 
veloped in  connection  with  soil  particle  acceleration  determination 
by  the  velocity  field  method.  This  concept  establishes  a relation 
between  the  shear  stresses  at  the  tire-soil  interface  and  the  mini- 
mum value  of  slip  that  is  necessary  to  develop  these  shear  stresses. 
Slip  that  exceeds  the  value  of  contact  slip  is  a waste  in  energy 
and  tire  wear. 

It  is  recommended  that  the  theory  of  contact  slip  be  further 
developed  and  general  relationships  between  contact  slip  and  inter- 
face friction  be  established.  It  is  also  recommended  that  in  all 
tire  tests  the  horizontal  displacement  of  selected  soil  particles 
at  the  surface  be  routinely  determined  and  compared  with  the  mea- 
sured conventional  slip  value.  The  horizontal  displacement  of  a 
soil  particle  at  the  surface  is  a measure  of  the  contact  slip.  Its 
measurement  is  a simple  way  to  obtain  experiirental  information  that 
could  lead  to  economies  in  off-road  vehicle  usage. 
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A towed  pneumatic  tire-soil  model  has  been  developed  that 
predicts  towed  force  coefficients  with  good  accuracy.  For  improve- 
ment of  the  model  it  is  recommended  that  tnree  dimensional  condi- 
tions be  taken  into  account  in  the  model.  The  consideration  of 
three  dimensional  conditions  is  also  essential  for  the  prediction 
of  towed  force  in  the  case  of  dual  tires. 
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APPENDIX  A 


COMPUTER  PROGRAM  FLOW  CHART 
FOR  ANALYSIS  OF  EFFECT  OF  SOIL  INERTIA  FORCES 
ON  TIRE  PERFORMANCE  BY  METHOD  OF  VELOCITY  FIELDS 
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APPENDIX  B 

COMPUTER  PROGRAM  FLOW  CHART 
FOR  COMPUTATION  OF  TIRE  PERFORMANCE 
WITH  CONSIDERATION  OF  INERTIA  FORCES 
BY  PARTICLE  PATH  METHOD 

(Pages  16,  17,  and  27  omitted) 
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0.  «UTOFLOlf  CHABT  SFI  - FlOCMAPT  ATIBE 


listing 


PAGE  OJ 


chart  TIT(_g  _ InTBOOUCTOFY  comments 


this  pbogpam  computes  off-boao  tiff  PFPFOBMANCe 

AT  VARIOUS  VELOCITIES  WIT„  CONSIOEPATI ON  OF 
SOIL  INERTIA  FORCES  COMPUTED  BY  THE  PARTICLE 
path  hETiiOO.  the  PROGRAM  CONSISTS  OF  THE 

follcking  parts 

MAIN  PROGRAM  IKTIREI 

subroutine  ISLFIC) 

SURBOUTINE  (KACCl 
subroutine  (PLOTSLI 
inputs  in  the  main  program  ABE  P£AO 
FROM  data  file  tyb.oat  In  the  follcming 
order 


CGPfCONE  index  GRA0tENT(L8S/CU  INI 
ci»coNf  Index  ilbs/so  fti 
RAsTiRf  Radius  ifti 
Qo=tirf  «idih  (fii 

LOtLOAO  (LBS) 

PO=L|M|T  PRESSURE  (LBSFSO  IN| 
OExOEFlECIIQn  coeff,  epsilon 
SLsSLIP  (decimal) 

SjsSLIP  parameter  J TERO 
Skksslip  parameter  k 
CFsCOHESION  in  front  field  (LBS/SD  FT) 
CflxCOHESION  IN  REAR  FIELD  (LBS/SO  FT) 
FFsFRIctION  angle  in  front  FIFLD  (DEGREES) 
FRsFRICTION  angle  in  bear  FIFLO  (DEGREES) 
GFXUNIT  meIGHT  in  front  FIELD  (LBS/CU  FT) 
Gfi.-UNIT  MEIGHT  IN  PEAR  FIELD  (L8S/CU  FT) 


IF  THE  PROGRAM  IS  USED  WITH  CONE  PENETEATION 
values.  CGR  or  Cl  IS  ENTERED  WITH  ThE 
other  as  TERO.  cohesion,  FRICTION  angle  and 
UNIT  WEIGHT  IS  COMPUTED  IN  THE  PROGPAM  AND 

Their  values  in  the  data  file  disregarded, 
alternatively,  if  The  program  is  used  with 

COHESION,  FRICTION  angle  AND  UNIT  WEIGHT 
VALUES  IN  THE  LAST  SIX  PLACES  OF  THE  DATA 
file  both  CGR  AND  Cl  MUST  BE  ZEPO, 
the  following  input  data  APE  ACCEPTED  IN 
The  RBOGRAM  from  TYPED  IN  DATA. 

VELwTRANSLATIONAL  VELOCITY  (FT/SEC) 

ZOERTM.OEPTH  of  influence  (FT)  —IF  NO  OTHER  DATA  ASSUME 

COITAL  TO  TJPF  PAOMJS 

zcta.factor  fop  particle  path  geometry— if  no  other  data 

ASSUME  AS  0,9 

ao-entry  angle  in  radians 


principal  VARIABLES  IN  THE  PROGRAM  ARE 


i 


j 
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AUTOFLOtf  CHADT  SET  - FLOCHAPT  KTIPE 


listing 


PACE  02 


CMAPT  title  - INTPODUCTOFY  COMMENTS 


AMkSEPAPATION  angle  in  paoians 
AIsTANCENT  of  log*  SPIAAL 
OtslNTFPFACE  FPtCTko.4  angle  IN  PADIANS 
LsLENGTH  OF  PASSIVE  /ONE  AT  SUPFACE 
SSL  IPs AVFPAGL  value  of  CONTACT  SLIP 

AZEPOsVfPTICLE  distance  BETvEEN  END  POINTS  OF  PAPTICLE  PATH 
CMEI«OME3»OME4sOMEGA  VALUES  AT  0,D«E  POINTS  OF  PAPTICLE 
TIMUTIMa.TlMAsELAPSEO  TIME  ASSOCIATED  FITH  BtO*£  POINTS* 
PESPECTI  VELV 

A( I • J)  «7( I *J>sX*7  COOPOINATES  OF  SLIP  LINE  FIELD  WITH 
PCFCPENce  TO  entpy  point 

S(l •J)«T(1*J>s$IGMA*  theta  value  at  I*J  point  Of  SLIP 
LINE  field 

A< Jt«0< J>*£( JIsCEnTPAL  ancle*  NQPMAL  STPESS*  SHEAR 
STRESS  AT  POINT  J OF  INTERFACE 
U(J)«V(J)sx  AND  2 coordinates  OF  INTERFACE  WITH  REFERENCE 
TO  T|P£  AXIS 

ACCX< 1 « Jl «ACC2( 1 • J >sX*7  COMPONENT  OF  ACCELERATION  AT  1,J 
POINT  OF  SLIP  LINE  FIELD 
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-}  I 
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/ PfAO  TQOM  oev  / 
/ 2 / 

/ VIA  format  / 
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/ INTO  TmE  list  / 

I 

I 

/XTI00320/I  NOTE  02 


• • 


• 4 4 • 


• LIST  4 CCP*  Ct« 

• PA,  BO.  LO.  PO. 

• OE.  SL.  SJ.  SKK, 

4 CF.  CPFF,  FP.  CF, 
4 OQ 

I 

I 

I 

/KTI00320/4  OJ 


• End  OF  DATA 


/ 7 / 

) I 

I NOTE  04 

44444444444 

• CONTINUE  4 

I 

I 

I 

COMPUTE  COHESION  AND 
fpiction  fopm  cone 
PENETPATION  DATA 
I 
I 
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CGP  .EO*  0 
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1 < 7 I « I 4AL0G10(CGP|  I 
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4 
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4 

4 

4 
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ITPUE 

1 

1 

! 

I /KTI 00*00/1 

O 1 

1 •--- 

PELO  4 1.1 


• I 
-M 
I 


FO  4 

57.3«ATAN< 1/<  U 
64  - .6e*PEL0)) 

GO  4 104  - 
17,54M  - PtLOl 

PELOP  4 PELO  4 
P0*< I - PELOI/SO 


/KT  I 004«.0/4 


PELOP  • 

• 1.1 


OT  « 


• false 


• I 

1 
I 

• I 

ITPUE  I 

I I 

1 I 

1 I 

I ♦- 

I A 

/KTI00440/I  II  I 

• * I 

I PELOP  4 l.l  II 

I 

1 I 

1 I 


I 


Cl  .EO*  0 


IFALSE 

I 

I 

I 

t 

I 

/KT I00400/I 

CO  S M.54CI 
CCP  « CO 
FO  « .2f*CI 
FOP  * FO 
CO  = too  ♦ .14CI 
CCP  » CO 

7 


I 


03.I0P— >• 

I 

FOP  K 

57«34ATAN|1/(U 
64  - •604PELDPM 

GOP  B 1C4  - 
I7.54I 1 - PELOP) 

7 

I 

t 

• • • 

• • 

• 3 • l6  • 

• • 

• ••  to 


I 

I 


CO  ■ CF 
COP  4 CP 
FO  s FF 
FOP  » FR 
CO  3 gf 


GOP 

3 GR 

03.1  34— 

1 

)1 

to 

1 t( 

FAR  s 

FO/57.3 

I 

I 

t 

/KT  130630/4 


FOP  #LT*  FO 


• false 


ITPUE 

I 

I 

I 

1 

1 

/KTI00630/I  la 



I FAP  3 F0P/57*J  I 

• , 

I 

I ( 

I 

COMPUTE  MAX  VALUE  OF 

I ntcPface  friction 

angle  DELTA 
I 
I 

/KTI30650/1  to 

• * 

OMAX  3 

. >04ATAN(SIN(FAP) 

) 

55  3 EKP(-(SL  4 
SJ)/SKK) 

£m3B  3 t - ES 

I 
I 
I 

COMPUTE  INTERFACE 
FRICTION  ANCLE  DELTA 
FP04  SLIP 

I 

I 

/KTI00690/I  20 


01  3 

ATAN(EK0B40MAX) 
00  3 57*3401 


1 

1 

/KTI00710/I 


PRINT  It 


21 

/ 


/KTt007f  /I 

• 4 4 4 4 4 

LIST 


4 4 4 4 


NOTE 

4 4 4 4 

CO*  FO*  GO 

4 4 4 4 4 

I 

I 

t 

I / 


22 
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LISTING 
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ms 


mm 


CHAPT  title  - PPOCEOUPES 


mm 


i 

0 

mm 


£ 


4 

d ft} 


e V 


( 

A 

Vt  K 


< 

m 9 


03*22-  — I • 
/KTI00730/I 

01 

/ 

/ 

/ 

/ 

/ PRINT 

12  / 

/ 

/ 

/ 

/ 

I 

1 

/KT10C730/I 

NOTE  02 

A A A A A A 

A A A A A 

A list  a RA 

• BO*  A 

A LO*  PO  • 

I 

1 

/KTI00760/I 

03 

/ 

/ 

/ 

/ 

/ PRINT 

14  / 

/ 

/ 

/ 

/ 

1 

I 

/XTIOOT-^O/I 

NOTE  04 

A A A A A A 

A A A A a 

A list  a OE*  00  A 

I 

/KTI 00760/1 

OS 

/ 

/ 

/ 

/ 

/ PRINT 

16  / 

/ 

/ 

/ 

/ 

1 

I 

/KTI00760/I 

NOTE  06 

A A A A A A 

A A A A A 

A list  a Cl 

A A A A A A 

* CGR  A 

I 

/XT! 00800/1 

07 

/ 

/ 

/ 

/ 

/ PRINT 

16  / 

/ 

/ 

/ 

/ 

1 

I 

/KTiooeoo/t 

NOTE  06 

A A A A A A 

A A A A A 

A list  a FOR 

• GOR*  A 

A COR 

A 

1 

/KTI00630/1 

09 

/ 

/ 

/ 

/ 

/ PRINT 

19  / 

/ 

/ 

/ 

/ 

i 

I 

/KTI00830/1  NOTE  10 

• •••••••••• 

A list  s SL«  SJt  A 

A SKK  A 


/KT100850/I 

II 

1 

A A 

A AAO  «LT* 

AAAM  a *001  4 

A A 

1 PI  a 3*14159 

I 

1 

I 

I 

1 VEL  a 0 

I 

I 

A A 

I 40IF  a .15 

I 

1 

ITRUE 

I 

1 

I 

I 

I 

1 PI  a 144AP0 

1 

1 

I AO  a 5 

I 

• wA 



I 

1 

I 

I 

1 

SET  minIMv;m  values 
FOP  COHESION  AND 

friction  ancle 
1 
I 
I 

/KTI00910/A  12 


IFALSE 


CO 

COP 


13 

>-A 

1 

I 

I 


30 


I 

-I  1 
I 


FO  .LE« 


I 

-M 

I 


FO  8 
ITER  8 


I 

•I  I 
1 


--04.14A — 

50 

I ARMAX  s 40 
I 

1 AR  a 10 

I 

I ARMIN  a 0 

I 

I ARE  « 70/57*3 
1 

1 PSt  a CO/TANOIFO) 

I 

I 

I 

ASSUNE  SEPARATION 

angle 

I 

1 

/KTIOIOSO/I  17 

1 ANC  a ATANMPI  A 1 
1 PSn  ATANO(OOI/Pt)  I 
1 1 
1 AAM  a AMC  I 

- -----  — ------------ 

I 

12*03—)  I 

51  1 


16 

—A 

1 


09.16— M 

52  I NOTE 


CONTINUE 


19 

i A 


1 

I 

I 

ASSJME  ENTRY  ANGLE 
I 
I 

/XTlOllOO/l 

t AAO  a AAM 
$.•••••••••• 

I 

06*06A*«>| 

I 

54  A 
A 


ADir 


21 

A 


false 


/KTIOtllO/1 
I AAO  • AAM 


22 
• 001  1 


I 

u — 

I 

COMPUTE  9EETL  FOR  LOG 
SPIRAL 

I 

1 

I 

I / 

/ s*ot 
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AUTOFLOw  CHAPT  SET 


FLOCHAPT  KTIRE 


LISTING 


PAGE  OS 


chart 


TITLE  - procedures 


0A.23— !♦ 
/KTtOl I3C/I 


RET (DC  *AAM/?«AAC  • 
AO/ST* 31 

P2  s Bl 

BE  « 

ALOGCCOSIAAC) 

/COS( AAM/2) ) 

RMAX  B 

9£/(AAM/2  - AAO) 

- ------------------- 

I 

I 

I 

/KTIOtlTC/*  e? 

• • 

false 


-•  bl  «GT.  pMAX 


ITPUf 

I 

I 

I 

I 


/•CTICU7C/I 


1 

Bl  B BMAX 

I 

— 

1 

OA 

1 

AAl  B APS INiei I 

1 

1 

FFl  e FO/57,3 

1 

1 

001  * 00/57*3 

I 

I 

GGO  s GO 

I 

1 

CCO  B CO 

1 

I 

I 

/KTI01230/1  QS 

------------------ -- • 

PSI  * CO/TANO<FO  I 

02  » 

UELT<001 «FFl  1 

Th2  « PI /2  ♦ 

• S*(02  ♦ ODD  • 

AAM  ♦ AA| 
--------------------- 

I 

I 

/KTI0l2eC/l  06 

- ----------------- — -• 

QUI  B OUP( OfFFI ) 

SIG  « PSI/OUl 

EPl  « 

EPOIFFl «TH2«C) 

SlOl  > SI G«EPI 

I 

I 

/KTI0I300/I  07 

• 

0U2  « 

OUA(001«FK1  ) 

0U3  « 

SIGI •0U2AC0SI00I I 
- CCI 

PRO  B - 

PA«EXP(-BI AAO/57 • 

3) 

I 

1 


im 

I 

CHECK  FOP  NO  failure 
CONDITIONS 

I 

I 

I 

/KTI0I340/F  06 

A • 

• A 

• • TRUE 

• 003  •GE*  PPL  •• — — 


A A 
A 

ifalse 


-0S*06a*.)a 

lie  i 


NOTE 

A A A A A A 

CONTINUE 


1 

1 

I 

COMPUTE  SEP*  ancle 
AND  BETA  IF  REAP 
FIELD  ONLY 

1 

I 

/KTI02I60/1  IS 


AAM 


1 *5*AR/57*3 


/KTI013Se/ 

I XXX  B 

1 OOF 


BETIDE* AAM/2* AAM* 
AO/S7*3> 

BE  > 

ALOGICOSC AAM) 
/cost AAM/2) ) 

A A 

I 
I 

/KTI02I90/I 


BMAX  s 
BE/l-AAM/2) 


I 
I 
I 

COMPUTE  FRONT  SLIP 
LINE  FIELD 

I 
t 

/KTtOl^dO/1  10 


1 

1 

1 

/KTI02200/A 


FALSE  • 
——A  Bl 


•GT*  BMAX 


......... 

.—  A 

I 

A A 

1 I 

H 

I 

A A 

6 1 

SLFI 

M 

I 

A 

« 1 

H 

1 

ITRue 

0 1 

H 

1 

1 

1 I 

H 

1 

1 

>-•>  A 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

/XTI02200/I 

IS 

KTI0I390/A  It 

1 

A- 

-A 

• 

A 

I 

I 

Bl  B BMAX 

I 

A 

A 

I 

A- 

• A 

A 

A 

TRUE 

I 

I 

A HMf 1) 

*E0«  0 

-A 

• 

A 

1 

1 

19 

A 

A 

I 

A- 

-A 

A 

A 

1 

I 

X9 

B PROASIN(AAM) 

1 

ifalse 

I 

I 

29 

« PROACOS(AAM) 

1 

I 

1 

N|  « 0 

1 

1 

1 

N2  « 37 

1 

KTIOIAIO/ 

12 

I 

1 

V9  * 0 

1 

1 00  B PPL  - OOF 

I 

I N|  B jjl 

I 

1 

I 

ADJUST  ENTRY  ANCLE 
I 
I 
1 

/KTIOIAAO/B  13 
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* I 

• I 
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A « • • • 
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• 
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I 
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I 
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PACE  07 
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Ans(oo)  - 

• 0*05*RR1.  • 1 

• * 1 
• • 1 


( ») 


• • • • 

, 7 

. 16 


I 

I 


NOTE 
CONTINUE 


1 

I 

1 

LlNir  PRESSURE  IS  NOT 
CONTROLLING 
1 
I 

/KTI01780/I  15 
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PPL 


OOF 


• • * • 


I 


06.06*— >I 

ao  I NOTE 


• • « « 


« AAO  .EO.  ARE 


7 

14 


/KTIOlf 40/1 
- ------------------- 

AMO  = AAO  - AA  4 


QM  s (PPL  * 
OQP 1/2 


AAO  3 AM0*Q0/0*4 


t 

I 

I 

/•CT101670/4 

« 

false  • 


-07*06*— J* 


• AAC  «EO«  • TRUE 
• (AAM  ♦ 0*001  1 

• • t 

• • t 


ifalse  * 


I 

I 

I 

STORE  DATA  FROM  FRONT 
F IElO 

I 

1 

/KT101810/I  17 


V9  s VV(JJl) 
s UU( JJl  1 
<9  3 VVK 
Z9  X UUK 
N2  s 37  - NX 


I 

I 

/KTI01860/I  NOTE  16 


• • 4 4 


aEClN  DO  LOOP 
110  N S N2*  36 


• 00  .LT*  0 
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1 

« 16  « 

-)  1 

* 4 

1 

1 
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I 

I 

19 
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1 

I 
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1 4 

• 4 

• 4 

1 

1 

1 I HMHCNl 

4 MM(37  - 

1 

4 

1 

1 
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N) 

I 

iTooe 

1 

/KTI01700/I 
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1 

I 

1 

1 1 000(N> 
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( 

I 

1 

I AAO  s 

AAO  A 

AAO  4 1 

I I 

N1 
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1 

1 

I 
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I 

I 1 

1 

I 

1 

I 1 EEECN) 

■ EE(37  . 
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1 

I 
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1 

/•CTI01670/I 

05  1 
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I 

1 1 VVV(N) 
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1 
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10 
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N> 

1 

I 

1 

4 

4 

I 

t 

* 

I 

FALSE  4 

4 

1 

I 

I 

1 

♦ -——4  AAO 

• LT. 

AAM  4 

1 /KT101910/1 

20 

4 06 

I 

1 4 

4 

1 4 

«4 

4 4 

1 

1 4 

4 

I 1 UJUfNI 

- UUC37  - 

1 

4 * 

1 

1 

4 4 

1 1 

NX 

1 

4 4 

FALSE  1 

1 

4 

1 I 

1 

4 OD  .CT,  0 

4 . 

1 

ITRUE 

I 1 THETA(N)  s 

1 

4 4 

A 

1 

1 

I I THETOT  - N) 

I 

4 4 

i 

1 

1 

I 

.. . . 

«4 

4 4 

4 

1 

1 

I 

1 

I 

1 

1 

I 

1 

1 

ITRUE 

I 

X 

I 

1 

I 

I 

I 

I /KT101710/1 

11 

I 110 

A 21 

1 

1 

I 4 

1 

4 4 

I 

1 

I AAO  > 

AAM  4 

0.001  1 

1 4 

4 

I 

1 

I 4 — 

I NO  A 

4 

( 

I 

1 

♦ end  of  00  4 

/KT101680/I 

07  1 

f... ........ 

- I I 

4 

LOOP  4 

*... 

...4  I 

1 

4 

4 

1 AOO  3 0.0> 

I I 

4 

12 

4 4 

4 

• --4  1 

4 » 

4 

I 

I 

4 

4 

IVES 

I 

1 

4 

4 TRUE 

I 

4 AAO  .CE*  ARE 


IFALSE 

I 

1 

1 

I 

I 

I 


•I* 

1 


13 

->4 


• 4*21 • 

• • 

* • • 54 


B-8 


X 1/0I/T4 

CMAOT  title  - PPOCEOUPES 


T 


*6^ 


^ * 


07.21  — )• 
/KTI0194C/I  NOTE 

4 BEGIN  00  LOOP 


01 

• • 


4 • 


112  N 

• 4 


Nl 


* ! 

HHIN) 

« 0 

QOCN) 

* 0 

EE<N) 

= 0 

VV(N) 

■ 0 

UU(N) 

■ 0 

TKT(N) 

= 0 

1 

I 

I 

:2  4 

03 

4 

4 

4 

4 

4 

END  Oi’ 

00 

LOOP 

4 

4 

4 

4 

4 

4 

IVES 

I 

1 

t 

I 

I 

/KTI02020/I 

1 40IF  s AAO  <-  AAM 

I 
I 

/<T(02030/( 

/ 


OS 

/ 


PFINT  99 


/KTIC2030/I  NOTE  06 

• U$T  « L«  AAO*  • 

• AOIF*  0QQ<N2),  4 

• HHM(N2»,  JJl  4 

44444444444 

r 

I 

/KTI02050/I  07 


I 


LFPONT  « L 


I 


/<TI02060/I 

/ 

/ 

/ PAINT  305 


I 

OPT  FOR  PLOTTING 
PROGRAM 

I 

I 

/KTI0210C/I 


oe 

/ 


09 


/ READ  FPOM  OEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 306  / 

/ INTO  THE  LIST  / 

t 

1 


AUrOFLOV  CHART  SET  - FLOCHART  KTIPE  LISTING 


- — 4 

/KTI02100/I 


NOTE 

4 4 4 

IGO 


I 

I 

/KTI021 10/4 


-06*1 14 — 11 
I 


14 

--4 

I 


1 

1 

1 

COMPUTE  INPUTS  FOR 
pear  slip  line  field 

I 

05.20—)! 

^ 1 15 

4-  — — — — — -«--4 

001  = « 00/57.3 
AAO  » AR/S7.3 
PRO  s - 

PA4FXP(-bl424AD/ 
S7.3) 

PAMIN  s 
Z9/COS(AAO) 

I 
I 
I 


/KTI02350/4 

4 

4 

false 


PRO  .GT. 
RAMIN 


ITRUE 

I 

1 

I 

I 

I 

/KTI02350/I 


1 

I 

/KTI02450/I 


/KTI02120/ 


13  I 
12  I 
I.  I 
fO  I 
II  I 


12 

--4 

H 

H 

H 


I 

I 

I 13 


It 


4- 

I 

RSO  » RAMIN 

• -4 

I 

1 

1 

to 

VVK  . V9  - X9 
UUK  « U9  • 29 

I 

t 

I 

JJl  4 0 

I 

1 

I 

PRt  • FOR/57. 3 

GGO  4 GOR 

I 

I 

/KTI02410/I 

19 

I 

CCO  4 COR 

I 

I 

AAM  s - AAM 

I 

I 

AAt  s ARSINIBI 1 

1 

1 

I 

I 

4- 

TM2  » PI/2  ♦ 
.54102  4 001)  - 
AAO  4 AAt 

1 

I 

I 

4 4 

1 

OU]  4 QUPIO.FFl ) 

-4 

I 

4 4 

I 

PSI  4 

1 

1 

IFALSE 

1 

1 

CC0/TAN(FF1) 

I 

I 

1 

1 

SIG  4 PSI/QUl 

I 

I 

1 

EPl  4 

I 

I 

1 

4- 

EPOCFFl.THC.O) 

1 

. B • 1 4 . 

120 
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AUTOFLOH  CHART  SET  - FLOCHAPT  KTIRE  LISTING 


PAGE  09 


08.20— 1« 

/KTI02A90/I 


I SIGl  s SIG*EP1 

1 0U2  « 

I OUAIOOl.FFI) 

1 

I 0U3  « 

I SIG140U2*COS(OOn 
I - CCI 


/KTI02520/* 


TRUE 

• 


ifalse 

1 

i 

1 

I 

I 

I 

check  for  no  failure 

CONDITIONS  IN  REAR 
I 


/KTI02540/* 


FALSE  • 

<-  — • QU3 


• GE*  PPL 


ITRUE 


1 

I 

I 

/KTI025A0/I 


1 XXX  3 3 

-09.024— 1 { 

• 17  1 

05 

1 OOF  > 1 

I 

I 

COMPUTE  REAR  SLIP 
LINE  FIELO 

I 

I 

/KTI02S70/I  06 

M ( H 

to  I SLFI  H 

I • I H 

>0  I H 

M I H 

• - 

I 

1 

I 

/KTI02580/*  or 

4 • 


XXX  •CE*  2 


IFALSE 

1 

I 

I 

I 

I 

------- 


• TRUE 


126 


• • TRUE 
• QQF  *L£«  PPL 

• ♦ I 

• • I 

• • I 


IFALSE  • 

1 

I 

I 

I 

I 

KTI 02620/1  NOTE 

BEGIN  00  LOOP 
127  N s 1.  36 


9 

19 


130 

09 


I 

I(- 

1 


10 


false  • 


00(N)  «GE. 
• PPL 


ITRUE 

1 

1 

1 

1 

1 

/KTI02630/I 


I 


OQ(N)  « PPL 


It 

I 


127 


I 

-U 

I 


12 


END  OF  00 
LOOP 


IVES 

I 

I 

1 

1 

I 

I 

« « « 

9,19l 

•••  130 


NO 


12S 


OOF  ,lt, 

4 <*954PPLI  « 


IFALSE 

I 

I 

I 

I 

I 

1 


9*  19* 
•••  130 


TRUE 


122  I NOTE  lA 

BEGIN  00  LOOP  4 
124  N 3 !•  36  4 


I 

M- 

I 


HHCNI  3 0 
00(N)  c 0 
EEINI  3 0 
UU(N)  3 0 
WIN)  8 0 


124 


16 


* 4 NO 

END  OF  DO 

• LOOP  4 


/KT102730/1 

- — ... 

VVK  m 0 

UUK  • 0 

VVA  8 0 

UUA  * 0 

JJl  3 0 

PPL  3 OOF 


/KTI02790/I 

AAM  3 • AAM 
AAF  3 0 
AX  3 0 
AY  3 0 


I 

I 

I 

• • • 

• • 

• 4.19. 

• 52 


130 

I 


09.384 — ) I 
I 

19 

1 lrear 

* t I 

I 

/KT102850/1 

20 

/ 

/ 

/ PRINT 

/ 

/ 

/ 

/ 

99  / 

/ 

/ 

I 

I 

/KTI02850/I  NOTE  21 

444444««,9, 

4 LIST  3 L.  AAO«  4 

4 0QQ(N2)«  HHH(N2).  4 
* JJI  4 

I 

I 

/KTI32860/I 

22 

/ 

/ 

/ PRINT 

/ 

/ 

/ 

/ 

307  / 

/ 

/ 

I 

I 

OPT  FOR  PLOTTING 
PROGRAM 

I 

/KT102890/X 

23 

/ READ  FROM  OEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 308  / 

/ INTO  THE  list  / 


I 

/KTI02890/I 

4 4 4 4 4* 

• UST  * 

4 4 4 4 4 4 

I 

I 

/KTI02910/4 


true  4 

•—’4  ICO  *E0« 

4 


NOTE  24 

4 4 4 4 4 
ICO  4 
4 4 4 4 4 


2S 

4 

4 

0 


IFALSE 

I 

1 

I 

I 

1 

/KTI02920/I 

4«***«— 

13  I 

tt  1 PLOTSL 

I • I 
10  I 
It  I 


. I 

, 

131 

I 

2; 

■S3 

3 N2 

- jji  ' 

N« 

3 N2 

- • 1 

/KTI02950/I  NOTE  28 

44*  *4**444« 

4 BEGIN  DO  LOOP  4 
4 160  N 3 N3*  N4  4 

I 

I 

I 

I / 


/10*01 
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CHART  title  - PflOCEOUFES 


I 


I 

I 

/KTI02990/I  0 

VVV(N)  » - VV(N  - 
N3  «>  1)  «■  VVA  » 

X9 

UUUINI  s UU(N  - 
N3  f U - UUA  ♦ 

Z9 

THETAINI  * Pt  - 
THETIN  - N3  ♦ II 

I 

I 

I 


COMPUTE  LOAOfORAKBAR 
AND  TORQUE  BY 
INTEGRATING  INTERFACE 

stresses 

I 
I 

/KTI03040/I 


Yl 

V2 

Y3 

YB 

NTHREE 


0 

0 

0 

0 

N3 


04 

mm0 

I 


/KTI 03090/1  NOTE  05 


BEGIN  00  LOOP 
190  N » NTHREE* 
36 


10*13— I 
I 


U2 


> (QOQIN  - 
1)  4 
2«000(N| ) 
/<3A(Q00<N)  ♦ 
oooiN  - im 

VI  * VVVIN  - I ) ♦ 
(VVVCNI  - VVVIN  - 
t 1 1 *U2 

I 

1 


HHHlNl  » - HM(N  - 

1 I 

I 

Ul  • UUU(N  -lit 

N3  4 1) 

1 1 

I 

(UUU(N)  - UUUIN  - 

I I 

1 

U)4U2 

EEE(N)  a - £E(N  - 

1 I 

1 

N3  4 U 

1 1 

I 

H2  a (UUUCN  - 

1 1 

I 

1 I - 

OOOIN)  a O0(N  - 

1 1 

I 

UUU(N) )/(VVV(N)  - 

4* 

N3  4 1 ) 

1 1 

-•  1 

I 

1 

VVVIN  -in 

— )« 

/KTI03120/I 


HI  a ATAN(H2> 


I 

I 

/KTI03I50/I 



A3  * Ul  - 29 

A4  a VI  > X9* 

AS  a SQRT(A3«42 

A44*2) 

U3  s •S«(QOO(N| 

QooiN  - in 

•-  — — 

1 

I 

/KTI03I90/I 


4 4 1 

4 4 1 

4' 

1 

U4  S «S4(EEE(N) 

— 4 I 

4 1 I 

• ABSIONF  - f 
4 LO)  - 0*14L0 

4 4 ! 

I 

EEE(N  - 1 )) 

1 I 

4 4 

END  OF  DO  4 I 

I 

I I 

4 4 

4 LOOP  4 1 

I 

US  a (VVV(N)  - 

1 1 

4 4 1 

I 

VVViN  - 1) >442 

1 I 

4 4 1 

1 

t I 

I |4) 

4 1 

1 

U6  a (UUU(N>  - 

I 1 

I 

lYES  I 

1 

UUUIN  - in  442 

I 1 

I 

1 1 

I 1 

I 

1 

I 

1 

I 

I 

I 1 

/KTI03220/1 

to  I 

I 

I 1 

4- 

--4  I 

294  4 16 

U7  a S0RTIU5 
U6I 


U74<UA*SIN(H1 > ♦ 
U3«C0S(H1 n 


U7*(U4«C0S(HU  • 
U3*S1K(H1 I) 


/KT1032S0/I 
Y6  * 

U7*U3A<-A3*SIN 
(MU  ♦ 
A44Cos(Hi n 

Y7  « 

U7«U4«(A3*C0S(H1 ] 
- A4*SIN(Min 


Yl 

Y2 


Yl 

Y2 


Y4 

Y5 


11 

•O* 

I 

I 

I 

1 

I 

1 

I 1 

I 

I 

I 

I 

t 


1 

I 

/KTI03290/I 

I Y3  » Y3 

1 

I Ye  « Ye 


Y6 

Y7 


12 

I 

I 

1 


I 

I ^ 


190 


I 


END  OF  DO 
LOOP 


/KT103320/I 

• 

I RN«  a B04YI 

I 

I pTF  = RC4Y2 

I 

I BTT  a B04YB 



1 
I 
I 


NO 

I 

I 

1 

• • • • 

• 10  • 
• 06  • 


/KTI03350/* 


IS 


false  • 


LO  *LT,  PNF 


TTFUE 


4 

• 

PNF 

4 

4 


ITRUF 

I 

I 

I 

I 

I 

/KTI03360/I  17 

I APX  a AP  I 

I 

I 

I 

4 IB 


FALSE 

4-4 

I 

I 

I 

• • • • 

• 11  • 
• 17  • 


/KTI03370/1 

I AMX  a - AAM 


19 

— 4 

I 


— )4 

360  I 20 


/ 

/ 

/ 

/ 

/ PRINT 

365  / 

/ 

/ 

/ 

/ 

I 

1 

/KT103010/I 

NOTE  21 

4 4 4 4 4 4 

4 4 4 4 4 

4 LIST  a 

BNF  • 

1 

+ 4.  t 4.  * 

1 

/KTI33830/1 

22 

/ 

/ 

/ 

/ PRINT 

370  / 

/ 

/ 

/ 

/ 

I 

I 

/KT133830/I 

NOTE  23 

4 4 4 4 4 4 

4 4 4 4 4 

4 LIST  a 

BTF  4 

1 

I 

/KTI33850/I 

24 

/ 

/ 

/ 

/ 

/ PRINT 

375  / 

/ 

/ 

/ 

/ 

1 

I 

/KT1338S0/I 

NOTE  25 

4 4 4 4 4 4 

4 4 4 4 4 

4 LIST  a 

BTT  4 

I 

I 

I 

COM'^UTE  sinkage 


I 

I 

/KTIOSeeO/l  26 

--------------------- 


I SNK  a UUU(N31  I 
I 1 
I NSN  a N3  I 

« — 


1 

I 


/KTI039C0/1  NOTE  27 

44444444444 

4 BEGIN  00  LOOP  4 
4 no  N a N3*  36  4 

4*44444444a 

I 

12*36 1 1 

I 

4 20 

4 4 

<.  4 

4 4 true 

4 UUU(N)  •GT*  4-4 

4 SNK  4 I 

4 4 I 

4 4 1 

4 *••• 

■FALSE  • 12  • 

I * OS  « 

1 *38S 

I 
1 
1 


Il2*06l 

* • • 360 
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• ( — 
s I 


01 

I 


I 

I 

I 

/KTI036IO/* 


F4LS6  • • 

• LO  .GT.  2*0NF  • 


ITPUE 

1 

1 

I 

I 

I 

/KTI03610/I 

• 

T F8  s 2 


1 

-U 

I 


AP  s AP  4 4AF0 


I AAM  8 - (AAM  - 
1 .074FBI 

1 

1 

1 

/KTI03640/4 

05 

4 

4 

« 

4 

false  4 

4 

4 AP  ♦GT* 

APMAX  4 

0 

4 

0 

• 

0 

0 

0 

ITPUE 

I 

1 

I 

I 

/KTI03640/I 

0 

AP  3 APMAX 


I 

-)  I 
I 


AX  s 0 

AY  3 0 


ie#i84— )4 
I 


1 1 • 23*— 

206  1 NOTE  08 

• CONTINUE  • 

I 

I 

I 

ADJUST  PEAR  AND 
separation  angle  TO 

NATCH  LOAD 

1 

! 

1 

/KTI03510/A  00 


17 


TRUE  • • 

• LO  *Gt.  UNF  • 


IFALSE 

1 

1 

I 

1 

1 

/KTI03520/1  10 


I 


F0  s 1 


I 


1 

1 

/XTI03530/* 

• 

« 

FALSC  4 


LO  «LT« 
•5«6NF 


ITFUE 


/KTI 03530/1 
I Fe  8 2 

♦ )i 


12 

t 


1 AP  3 AP  - 4«F6  I 

I I 

1 AAM  s - <AAM  ♦ I 

I •074F6I  I 

-* 

I 

I 

I 

/KTI03560/4  14 


false  • • 

------  AP  *LT,  ARMIN  • 


ITPUE 

I 

I 

1 

I 

1 


/KTI03560/1 


false  • 4 

4— “4  LO  •GT*  6NF 


ITPUE 

I 

I 

I 

I 

I 

/KTt03360/I  16 

• - — « 

1 APY  b AF  I 

I 

) I 

I 

4 10 


false  4 4 

* LO  .CT*  BNF 


/KTI03300/I 



t AMY  s “ AAM 


20 

1 


1 

•)  I 
I 


4 AP  »C0*  APMAX  4- 

4 4 


false 


/KTIC3410/ 


4 AP  »E0»  APMIN  4- 


/KTI03420/4 


TRUE  4 

4 APX  *E0* 


I 

1 

i 

AP  s 

APMIN 

I 

IFALSE 

1 1 

AR  3 *54(ARX  4 

1 

1 

I 

- — 4 

1 I 

ARY) 

1 

I 

I 

I 

I 1 

I 

• * * 

4-- 

- --- 

- ----- 

1 

1 I 

AAM  s *54(AMX  4 

I 

• • 

1 

16 

I 1 

AMY) 

1 

• 12*01 . 

' •— 

— — — • 

1 I 

I 

• • 

I 

AX 

* 0 

I 

------------- 

— 4 1 

> 

M 

M 

O 

I 

• • • 300 

1 

1 

1 

I 

I 

AY 

B 0 

I 

1^ 

AY  ■ 0 

1 

1 

I 

1 

I 

330  I 24 

/ . 
/ / 

/ PRINT  340  / 

/ / 

/ / 

I 

I 

1 

• • « 

• * 

• 5*21* 

» • 

•••  909 


- — )• 

349  I 25 

/ / 
/ / 

/ PRINT  350  / 

/ / 

/ / 

I 

I 

I 

• • • 

* • 

* 5.2U 

• • 

•••  999 


•— 14 

I 

/KT103430/4 


ARY  .EO* 


/KTIG3440/ 


false 


TRUE 

4-4 

I 

I 

1 

• • • • 

» 11  • 

• 08  • 

*296 

27 


300 


300 
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chart  title  - PROCEOUPES 


9 *' 


1 1 .07* — II 

I 


I ITEP  s ITER  ♦ I 1 

*-->> — • 

I 

I 

I 

/KTI036R0/*  02 


• ITER  .GE*  20  •- 


/KTI03700/*  03 


• • TRUE 
XXX  .EQ*  2 

* • 1 

• • I 

• • I 

IFALSE  /’s*. 
I • 14  • 


329  I 04 

/ / 
/ / 

/ PRINT  330  / 


* 5«21* 

999 


/ 365  / 


10*26 M 

I 

I SNK  s UUU(N) 
I 

I NSN  s N 


10.26— > I 
I 

360  4 06 


END  OF  00 
LOOP 


lYtS 

I 

I 

1 

I 

I 

/KTI03960/I 


I SNKIN  a 124SNK 


NO 

I 

I 

I 

*10*. 

26  . 


07 

I 


I 

I 

/XT I03970/I 
/ 


06 

/ 


PRINT  390 


I 

I 

/XTI03970/I  NCTF  09 

• •«4**4***» 

• LIST  = SNXIN  • 

1 

I 

I 

COMPUTE  INPUTS  FOP 
particle  PATH 
GEOMETRY 

I 

I 

/XTI04000/I  IC 


PUTO  « UUUfN3) 


- VVVINSN)  I 
I 

• VVVIN3I  I 


I 
I 

I XOIS  4 
I 

I XPPOJ 


I 

I 

COMPUTE  PULL  AND 
TOPOUF  COEFFICIENTS 
I 
I 

/XTI04040/I  I 


PUN  « 0TF/PNF 


TON  * 

BTT/’|5&46NF) 


I 

I 

/XT  lOA 060/1 
/ 


12 

/ 


PRINT  39** 


I 

I 

/XTI04060/I  NOTE  13 

• *444444444 

• LIST  ■ PUN.  TON  • 

1 

I 


— ) • 

/KTI04060/I  NOTE  14 

• BEGIN  DO  LOOP  * 

* 444  N s N3.  36  4 


1 

1 

15 

/ 

/ 

/ 

/ 

/ PRINT 

433  / 

/ 

/ 

/ 

/ 

I 

1 

/KTI04090/I 

NOTE  16 

4 4 * 4 4 4 

* * * * 4 

* LIST  s N. 

HHH(N>.  4 

* JJO(N). 

EEE(N).  4 

* VVV(N). 

****** 

UUU(N)  4 

4 4 4 *4 

1 

I 

444  4 

17 

* 

* 

4 

* 

• 

o 

z 

4 

4 LOOP  * 

* 

4 

4 

4 

* 

I 

YES 

I 

1 

/XTI04120/1 

16 

/ 

/ 

/ 

/ 

/ PRINT 

309  / 

/ 

/ 

/ 

/ 

1 

1 

/XTI34120/I 

NOTE  19 

4 4 4 * 4 * 

4 4 4 4 4 

* LIST  » 

RUTO*  4 

4 XOIS* 

4*4*44 

XPPOJ  4 

4 4 4 *4 

1 

I 

/XTI34140/I 

20 

/ 

/ 

/ 

/ 

/ PRINT 

402  / 

/ 

/ 

/ 

/ 

I 

I 

I 

FOR  THE  COMPUTATION 
OF  INERTIA  FORCES 


I 

I 

/KTI04170/I  21 


/ READ  FROM  OEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 400  / 

/ INTO  THE  LIST  / 


I 

1 

/KTI04170/I  NOTE  22 

4*44*444444 

• LIST  S VEL*  ZDEPT  « 

4444**44444 

I 

I 

I 

I / 
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AUTOfLOW  CHART  SET  - FUOCHA'  ' KT|PE 
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chart  TITLE  - PROCEDURES 


/KTI04IR0/*  01 


VEL  .EO« 


TRUE 

1 

I 

1 


t^ALSS  • £»  • 

I • 21  * 

I • • • • 

I 999 


I 

/KTI04200/I  *^076 

44444«4«« 


4 

BEGIN  DO  LOOP 

4 

4 

4 

410  N a N3f  36 

4 

» 4 

.. 

I 

1 

03 

1 

tmetacni  * 

1 

«540I  - THETA(N) 

1 

XA(N)  s VVV(N)  - 

I 

X9 

1 

ZAfN)  a UUUINI  - 

1 

4- 

Z9 

1 

END  OF  00 
LOOP 


4 

IVES 

1 

I 

I 

I 

/XTI04260/I 

05 

I ChO  * 0 

I 

1 SSL  « 0 

1 

I CHOSL  a 0 

I 

I 

I 

/KTI04290/I 

06 

/ 

/ 

/ print  415 

/ 

/ 

/ 

/ 

/ 

/ 

/ 

1 

I 

/KTI043I0/I 

07 

1 N M J6  1 

I 

t 

1 

COMPUTE  VELOCITY 
BOUNDARY  CONDITIONS 
AND  AVERAGE  CONTACT 
SLIP 

I 

i3«ta— -n 
Ai7  I oe 

... — 

I RAO  ■ 

I S0RT(ZA(N)442  ♦ 

I XA(N)442) 

I AI4  > 

I ATAN(XA(N)/ZA(N|) 

I 
I 


/KTI04350/4 


«LT«  36 


IFALSE 

I 

1 

I 

I 

I 

/KTI 04360/1 


420 

I 

-M 

I 


A15  « - 
ATAN(CZA(N)  - 
ZA(N  ♦ 

tn/<XA(N)  • 
XACN  4 II ) I 


11 

I 

I 

I 

I 

I 


A15 


A|4 


A1 


10 

I 


430 

AKl 


I 

-)  I 
1 


c C0SIAI4  - 
A15I 

AKS  ■ 

TANITHETA(N) ) 
AK2  s 

AKS4SIN(AI5)  ♦ 

cost AtS) 

AK3  ■ 

VEL/I AK14AK2I 


12 

t 

I 

f 

I 

1 

I 

I 

1 

I 

I 

I 

I 


/KTI04430/1 

4 

! 


AI6  « l.O  - 
(Ak:  .'‘AK21 


I 

I 

I 

/KTt04440/4 


•EO*  36 


iealse 

I 

t 

I 

I 

t 

/KTI044SC/I 


IS 


CHORD  * 
SORTICXAIN)  • 
XA(N  - 1 M442  ♦ 
(ZA(N)  • ZAIN  - 
1 ))442) 

CHOSL  « AI64CH0PD 
SSL  » SSL  ♦ CHOSL 


I 

I 

/XT104480/I 
I CMO  ■ CMO 


16 

-------- 

CHORD  I 


I 

II- 

I 


N - 1 


17 

I 


I 

I 

I 

/KTI04510/4 


IB 


N .CEt  N3 


IFALSE 

I 

I 

I 

I 

I 


TRUE 

4-4 

I 

I 

I 

• • • • 

• 13  • 

• OS  • 


— ————14 

/KTI04S20/I  19 

•—•——•4 
1 3SL1P  « SSL/CHO  I 


I 

I 

/KTI04S30/1 

/ 


/ PRINT  460 


I 
I 

/KT104530/I  NOTE  21 


4 4 4 4 


4 4 4 4 


4 LIST  s LFRONTf 
• LREARf  SSLtP 

I 

I 

/XTI34550/1 

/ 


22 

/ 


/ PRINT  465 


I 

I 

I 

ASSUME  ZETA  FOR 
PARTICLE  PATHIAS  0*9 
IF  NO  DATA) 

I 

1 

/KTI04580/I 

/ READ  FROM  DEV  / 
/ 5 / 

/ VIA  FORMAT  / 

/ 470  / 

/ INTO  THE  LIST  / 

I 
( 

/KT104580/I 

4 4 4 4 4 4 
4 LIST  « 

4 4 4 4 4 4 

I 
I 
I 

COMPUTE  RHO  FOR 
PARTICLE  PATH  FROM 
SLIP 

I 

/KTU4610/I  25 


NOTE  24 

4 4 4 4 4 
ZETA  4 

4 4 4 4 4 


ZKSQ  > ZETA442 
ZKCU  • ZETA443 
ZXFO  ■ ZETA444 


I 

I 

/KT104640/I 


ABO  • - 

t 

liAOl l**ZKSO 

•4 

1 

2*372184ZETA 

4 

1 

2«0t924 

1 

ABl  • - 

1 

•2737144ZKS0 

4 

X 

U050394ZETA 

t 

4- 

2«0187 

X 

I 

I 

/KTI04660/I 


I 

A82  ■ - 

-4 

t 

10.9S64ZKCU  4 

X 

24.66444ZKS0  - 

17.62944ZETA  4 

X 

4- 

A.71601 

I 

I 

I 

1 / 
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I 


/ 464 


/ 


13*?T- — I*  ♦ 

/•CTI04670/X  01  I 

• - 1 

I Ab3  s > 11 

1 717.661*;kFO  ♦ I 1 

I 250«*05«ZKCU  - I X 

X 3263*0«2KSO  ♦ 1 I 

I X663*36«2ETA  - 1 I 

1 396.737  I 1 

• t 

I X 

I 1 

/KT1046P0/I  02  I 


1 QHO  r 

X AB3*(SSLIP««3I  * 
I AB2*1  SSL IPA*2)  ♦ 
I AR|*SSLtP  ♦ AOO 


I 

t 

/ICTI04700/X  03 


/ / 
/ ✓ 

/ PRINT  100  / 

/ / 

/ / 


1 

X 

/KT104700/I  NOTE  04 
*••••• 

• LIST  « ZCTA.  • 

♦ SSLIP.  PMO  • 

I 

I 

/KT104720/I  06 

--------------------- 

OM7  s PH047ET4 

FSTAO  = 

ATANl .642ETA/TAN 
(.254P1  - 
.54FF1 ) I 


I 

I 

I 

COMPUTE  OMEGA 
PARAMETER  ANO  TIME 
VALUES  for  particle 
PATH 


/KTI04750/I 


06 


0M£0  a 

2«FSTAR/(t  > QHOI 

OMFl  a Pl/l 1 ♦ 
PH042ETA1 

OME2  « 24PI/U  ♦ 
«HO) 

OMES  « 34PI/U  ♦ 
OH042ETA ) 

I 
X 

/KTI04790/1 

OMEVE  a 
24P|/<PH0  - 1 I 

OMTEN  a 

• 5«(0MEVE  f 031 

QEVE  a OMEVE 

TIMO  « 0 

003  • 03 


I 

I 

/KT104640/I  06 

• c... 

I lOUlT  • 0 1 

• 

I 

I 


— — — 

472  1 

09 

X 

a- — 

. _ _ a 

X 

1 OFNO  s 

1 

1 

1 .SaZETAasiOlPHOa  1 

1 

I /FTA.OMTEN) 

t 

1 

X /COC(PHOtOMTEN)  1 

X 

•—a 

1 

X 

I 

1 

X 

1 

X 

/xTX04660/a  to 

X 

a a 

X 

a a 

X 

aABSIOENO  - a 
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X 

a 1.0)  .LT.  1.0 

a 

--- 1 - 

a a 

X 

474  1 

16 

a a 

1 

a 

a a 

1 

X 0MF4  « OMTEN 

X 

a 

X 

I 

I 

ifalse 

I 

1 TXMl  r LFPCNT/VEt  X 

I 

1 

1 T|M3  * TXMl  ♦ 

1 

I 

I 

X XOIS/VEL 

X 

1 

1 

1 TXU4  B T|M3  f 

1 

1 

X 

1 (XPROJ  - XOXS 

/KT|04870/a  tl 

t 

1 LRFAP)/VEL 

X 

a a 

1 

a-o> — — .... 

— a 

a a 

X 

I 

false  a a 

X 

X 

a 

X 

/XT I0499C/X 

19 

a a 

1 

a^'o  — — — — — O'-— — • 

--a 

a a 

I 

I A2EPO  : 

I 

a a 

1 

1 PUTD/(COFtPH2» 

X 

a 

I 

X OMEO)  - 

I 

XTPUE 

X 

1 C0E(PH2»0MFA)) 

X 

1 

1 

— a 

X 

1 

I 

1 

1 

-------- 

/XTX04670/1 

12 

X 

« 

— * 

t 

I 003  a OMTEN 

I 

— - • 

X 

1 

a 13 

X 

a a 

1 

a a 

X 

FALSE  a a 

X 

a 

X 

a a 

X 

a a 

I 

a a 

X 

a 

X 

XTPUE 

I 

I 

1 

X 

X 

1 

/XTX  04680/1 

(4 

X 

I 

X 

1 QEVE  s OMTEN 

t 

X 

— --- — 

1 

--a 

I 

1 

1 

15 

1 

X 

— ....... 

--a 

X 

X lOUfT  S XOUIT  ♦ 

I I 

X 

— ....... 

--a 

X 

X 

/XTX  04900/1 

16 

1 

X 

X 

X 

X 

/KTX04920/a  20 

/ 

/ 

1 

X 

a a 

/ 

/ 

1 

1 

a a 

/ PRINT  476 

/ 

1 

1 

a a 

TRUE 

/ / 

1 

I 

a XOUXT  .GT.  20 

a-a 

/ / 

1 

X 

a a 

X 

1 

1 

X 

a a 

I 

I 

1 

X 

a a 

I 

1 

1 

1 

a 

/KTl 04900/1  NOTE 

17 

1 

X 

XFALSE  • 

5 • 

aaaaaaaaa 

a a 

X 

X 

X • 

21  a 

a list  « Umten. 

a 

I 

X 

1 

a OENO 

a 

a a 

1 

1 

X 

X 

1 

1 

999 

! 

1 

1 

X 

X 

X 

/XTX04930/I 

21 

I 

X OMTEN  B .54(003 

♦ 1 

X 

X OEVE) 

I 

I 

1 

a— 

X 

--a 

X 

X 

1 

-)  I 
1 


/ 

/ 

/ PRINT 

/ 


460 


/ 

/ 

/ 


/ 


/ 


I 

I 

/KTI06910/I  NOTE  23 

• aaaaaaaaaa 

♦ LIST  s OMEO.  a 
4 OMSt • 0H£2«  0ME3.  * 

* 0ME4  a 

X 

I 

/KT105030/I  24 


f 

/ 

/ PRINT 

/ 

/ 


461 


/ 

/ 

/ 


/ 


/ 


1 

1 

/KT10S030/I  NOTE  25 

aa«  44A4«««« 

a LIST  s TXMl • a 
a TIM3.  TIM4  a 

I 

1 

I 

/KTI05050/a  26 

a a 
a a 

a a TRUE 

a <XX  .EO.  2.0  a-a 

a a 1 

a a 1 

a a 1 


FALSE 


• 15  • 

• 03  • 


/KTI05060/ 


CXI  a OMEVe/PX 

RAT4  « 

tanccohes  - 

•5aoM£VE)/C<t) 

RATI  a 

TAN( (OMEO  - 

•SaOMEVE )/CKl) 

I 

I 

/XTI05090/1  26 

a ----------..--------a 

I PAT2  - 1 

1 RATA/QATl  - 1.0  I 
I 1 

1 TIMO  a TIM3/RAT2  1 
I 1 

I CX2  * 1 

I RAT4/<T1M3  ♦ 1 

1 TIMO)  I 

1 1 
I OMCt  a I 

I CKtaATAN(CK2a  1 
1 TIMO)  a .SaOMSVE  I 


I 


1 

/KTUS130/I  29 

a— 

1 0MC3  • X 

I CKL«ATAN(CK2a  ) 

1 (7IM3  a TIMO))  a I 

1 .SaOMEVE  1 

a---- 


X 

X 

X 

X / 


/ts.ot 
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AUTOfLOW  CMAPT  SCT  - FLOCHAPT  KTIPF 


LISTING 


PACE  15 


CHAQT  title  - PROCeOUPCS 


14*29  — M 

/KTI0M40/I  01 


/ / 
/ / 

/ PRINT  4FA  / 

/ / 


/ 


/ 


I 

I 

/KTI05140/I  NOTE  02 

*44«»4«4*4« 

• LIST  a A7EP0,  • 
4 TIMO*  T1M3*  TIHP,  4 
4 OMCl*  0MC3  4 

I 

14*26—1  I 

402  I NOTE  03 

44444444444 

4 CONTINUE  4 


I 

I 


COMPUTE  CONSTANTS  IN 
OMEGA-TIME  RELATION 

I 

/KT I 05 160/ I 

04 

1 CKl  4 OMEVE/PI 

I 

1 PTA  a TAN( (0ME4 
1 •S4CWEVE)/CK1 1 

1 

I PTI  a TANUOMEt 
I •S4CMEVE)/CK1 ) 

1 

I 

I 

/KTI05210/I 

05 

1 PT2  a PT4/PT1  - 
1 1 *0 

I 

I 

r T|MO  a TIMA/PT2 

1 

I CF2  a PT4/<T|M3 

1 TIMO) 

4 1 

I 

I OMCl  a 

1 CKI4ATAN(CK24 

1 TIMOI  4 *540MCVE 

I 

I 

I 

I 

I 

/KTI05290/I 

06 

I 0MC3  4 I 
1 CK14ATAN(CK24  1 
1 <TIM3  ♦ TIMOI I ♦ I 
1 «940MEvE  I 
1 1 

AZEPO  a 

PUTO/<COE<PHZ, 
OMEU  - 
C0E<PH/«nME4U 


I 


I 


/KT 109270/ I 

07 

/ 

/ 

/ 

/ 

/ PRINT 

466  / 

/ 

/ 

/ 

/ 

♦ 


I 

/KTI09270/I  NOTE  06 
44444444444 

4 LIST  4 A7EP0*  4 
4 TlMCc  T|M3,  TIMR*  4 
4 OMCl*  0MC3  4 


I 

I 


4 

/KTI05260/I  09 


4--- 

AX  a 0 

— •4 

AY  • 0 

ARX  a 0 

AMX  a 0 

APY  a 0 

4--- 

AMY  a 0 

• — 4 

I 

I 

I 

RECOMPUTE  SLIP  LINE 
FIELDS  FOR  INERTIA 
FORCES 

I 

I 

t 10 


• A • 1 6* 


• • 

50 
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AUrOFLOV  CMAOT  SET  - FLOCHAPT  KTIPE 


LISTING 


PACE  11 


CMAPT  TITLE  - SUHPOUTINE  SLF I 


I 

-U 

I 


NOTE  01 


BEGIN  00  LOOP 
3 I « 1 « 49 


4 4 4 


I 

n 

I NOTE 

4 4 4 4 4 4 4 

BEGIN  00  LOOP 
4 J a U 17 


02 


I 

->  I 


X(  V • J > a c 
7( I tJ)  a 0 
S(  I .Jl  a 0 
T(t,J)  = 0 
I 
I 


04 


FNO  OF  no 

LOOP 


ENO  OF  00 
LOOP 


lYES 

1 

I 

I 


I 


/SLF002S0/I  NOTE  06 


BEGIN  00  LOOP 

6 J a } , 36 


1B«09*- 


I 

•II 

I 


A(  Jl 
B(  Jl 

C(  J) 
0<  Jl 
U(  Jl 
V(  Jl 


I 


/SLF00320/I 

4 — 

H(  J) 
0(  Jl 
E(  Jl 


->4 

c 


ENO  OF  00 
LOOP 


/SLF00420/I 

OAF  s TANCOll 
II  « 24N  - 1 
Cl  a CO/TF 
Ol  a PI/4  - Fl/2 
F3  a t - SINIFI I 


/SLF004  70/I 

V3  a - TANIFl  > 

V4  a 1 

VS  a TANCFll 
V6  a 1 
02  a 0 


/SLF00520/I  13 

--------------------a 

03  a OEl(02*F1) 

T|  a (03  - 021/2 

OS  a 014( 1 4 
OFI/2 


04  a 0EL(0S*F1 I 


/SLF 00560/ I 


lA 


• 


T2  a p(/2 
• S«(04  ♦ 01  I - 
AO  ♦ A1 


I 
1 
I 
I 
I 
I 

T3  a •7S»P|  4 I 
•54FI  - T2  4 T1  • I 
AO  4 Al  I 


VI  a Cl 


) a 

/SLF00590/I 


IS 


VEL  .EO* 


TPUE 


/SLF00660/ 


19 


LI  a 11 

DISC A0*AM,T2»T1 • I I 

FI  I II 

I 1 

OMAX  a I ( 

ATAN(SIN(F1I)  I 1 

I I 

FAC  a 2 4 I I 

•5A01/0MAX  I 1 

------- — — ...«  J 

I I 

I I 

/SLF00620/I  16  I 

•  4 I 

I L2  s 11 

I FAC4L14C05C U1 ) 1 I 

I 4C0S(T3  - I I 

I Ft)/C0S(F1I  I I 

• 1 

I I 

I 1 

1 I 

ESTIMATE  L DISTANCE  I 

FOP  PASSIVE  ZONE  1 

I 1 

I 1 

/SLF00640/I  17  I 

• - — 4 I 

I L a - t2  I I 

• - • I 

I I 

I I 

I I 

/SLF006S0/4  le  I 

a 4 I 


XX  ,E0«  UO 


ITFUE 

I 

I 

I 

I 

1 

/SLF00670/I 


I 


22 

— 4 

1 


•19*164- 

20 


23 


/SLF00690/I 

/ 


24 

/ 


/ / 
/ PRINT  19  / 

/ / 

/ / 


/SLF39690/1  NOTE  25 

44444444444 

A LIST  3 L«  AO  4 


SET  dOUNOARV 
CONJITIONS  AT  FREE 

surface 

I 

I 

/SLF00720/I  NOTE  26 
44444444444 

4 9P.GIN  00  LOOP  4 
40  J a 1«  N 


4 4 4 


19*09- 


I 

-II 

I 


I a N 4 1 - J 

AJ  a J 
AN  a N 


Z(  I *J> 


X(  I *J|  a (AJ  - 
2>4L/(AN  - 2) 


/SLF30760/I  29 



1 X(N»1I  a 0 t 

I I 

I X(N  • 1 *21  a I 

I *54L/(AN  - 21  I 

* — 

I 

I 

I 

/SLF30900/4  29 


4 VEL  *E0*  0 


TRUE 

4-4 

I 

1 


IFALSE  • 19  * 
I * 07  • 

1 • * • • 

1 22 
1 
1 
I 

COMPUTE  INERTIAL 
ACCELERATIONS  AT  FREE 
SURFACE 

I 

I 

/SLF00820/I  30 


12  1 
I 6 I 

I • 1 
10  1 

II  1 


- -4 


ACCE 


L a lPEAP 

I 

I 


/19.01 


X x/o\/y* 


AUTOFLOV  CH>PT  SET 


FLOCHAPT  KTIPF 


listing 


PACE  19 


CHART  title  - Subroutine  slf r 


IB.31— -)  • 

/SLF00630/t  01 

« 

02  s 1 

ATANIACCXI 1 • J)  I 

/IACC7(I»J>  *■  I 

COM  1 


I 

I 

1 

/suFooeAo/* 


AeSIOSI  *GF« 
•A6S< •99*rt !• 


iTPur 

I 

I 

I 

I 


SLf  OOBAO/t 


02  8 

«og«Fl*02/AeS(02l 


03 

- • 

I 


t 

O I 
I 


ifalse 

I 

I 

I 

I 

I 

/SLF00660/t 


06 

/ 


PRINT  360 


/SLf00e60/l  NOTE 
• L I ST  « j«  02 


06 


--I  e*29«— » I 
22  I 


03  5 0FL<02*F| ) 
Tl  » (03  - 021/2 
lr2  s QUR(02«F|  | 

SI  s «l/«2 

sn.ji  s SI 


I 

I 

/SLF 00930/1 



I TU»JI 


oe 
--  • 
t 


I 


FND  OF  DO 
LOOP 


IVES 

1 

I 

I 

I 


NO 

t 

1 

I 

« • • • 

• 16  • 
• 27  • 


/SLF00950/I 

• -- 

t NET  = 2*N  - 


1C 

I 


1 

t 

1 

COMPUTE  SIGMA  AND 
THETA  AT  SINGULAR 
POINT 

I 

I 

SO  I NOTE  It 

• PEGIN  DO  LOOP  • 

• 70  I 8 N*  NFT  • 


I 

•)  I 
I 


J s I 

At  = FLOAT!  I I 

AN  s FLOAT(N> 

X( !• J)  « 0 

7( l»J)  * 0 

I 

I 

/SLFOI030/I  I 

Tl  1*3)  **  Tl  ♦ 

(72  • Tll*|A|  - 
AN)/(  AN  - I ) 

POM  • 

rPO(F|*T(|»J>*T|) 

S(  I • J)  • 
P0W»S(N*1  ) 

I 

I 


END  OF  DO 
LOOP 


I 

I 

1 

1 

/SLF01070/I  i 

• 

I 02  8 OUAlDt  »FU 

Oil)  » 02ASI 2AN  - 
t*l)*COS(Dn  - Cl 

Ell)  8 (0(1)  * 

Cl )*DAF 


/20.0I 


ll/Ot/74 


AUTOfLOW  CH»PT  SET  - ELOCHAPT  KTIPC 


LISTING 


PAGE  20 


Si 


CHADT  title  - SU0POUTINE  SLF I 


10.16-— )• 

/SLECI ItO/I 


I THETU)  » T<2*N  - 

1 1 .n 

I 

I J S 1 


COMPUTE  COOPOINATES* 
sigma  and  theta 
VALUES  FOP  SLIP  LINE 
FIELD 


I 


/SLFOtlAO/l  NOTE  02 


BEGIN  00  LOOP 
200  J s 2*  N 


I 

— M 
I 


JNT  » J ♦ 2*<N  - 
1 I 


26.12*- 
60 


I 

>)  I 
I 


BEGIN  00  LOOP 
100  I s KNT«  JNT 


•H 

I 


06 


1 •t0«  <J  ♦ 

•2*<N  - in  * 


TPUE 

I 

I 


/SLF0I20C/ 

K 


FALSE 
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